THE BULLETIN OF 


ae 


Mathematical 


BIOPHYSICS 


DECEMBER 1959 


A Mathematical Interpretation of Motor Reactions to Fixed and Mov- 
ing Stimulation of the Eye—Maximo Valentinuzzi - - - - - - - 327 


A Model for Unequal Ventilation of the Lungs Assuming a Common 
Dead Space and Two Separate Dead ee E. Wise and J. 


G. Defares - - = --- +--+ 2-22 2s ee eee tee 343 
: Stochastic Models for Chemical Reactions: II. The Unimolecular 
Rate Constant—Anthony F. Bartholomay- - -----+-+-- - 363 
Mathematical Biophysics of Automobile Driving—N. Rashevsky - - 375 
Ona Logical Paradox Implicit in the Notion of a Self-Reproducing 
: Automaton——Robert ebe e ale ins fle) 6 a0 me ie = 6 eis 387 
| "Mathematical Biophysics of Color Vision: II]. Color Constancy— 
395 


_ WH. D. Kandohl- = 3 = = + +--+ 6 Se ee ee 


ng is +1960 International Congress for Logic, Methodology and Philosophy 
a te iba Science: [RO ees ea ae ee ee a 


- 
—_ 


- CHICAGO 
NUMBER 4 


The Bulletin of 
MATHEMATICAL BIOPHYSICS 


Eprror: 


N. RASHEVSKY, UNIVERSITY OF CHICAGO 
AssociaTE EpitTors: 
H. D. LANDAHL, UNIVERSITY OF CHICAGO 
ANATOL RAPOPORT, UNIVERSITY OF MICHIGAN 


The BULLETIN is devoted to publications of research in Mathematical 
Biology, as described on the inside back cover. 


AUTHORS: 


Please read carefully the inside back cover. 


Published March, June, September, and December. Subscription 
rates per volume: U.S.A., $10.00; Canada and Pan America, $10.50; 
all other countries in Postal Union, $11.00. Single copies: $3.25. 
Checks should be made payable to The University of Chicago Press, 
5750 Ellis Avenue, Chicago 37, Illinois. Subscriptions will be en- 
tered to start with the first issue to be published after receipt of 
order. Editorial correspondence: The Editor, The Bulletin of Math- 
ematical Biophysics, 5741 Drexel Avenue, Chicago 87, Illinois. 
Subscribers are requested to notify the Press and their local post- 
master immediately of change of address, giving both the old and 
the new address. 


Second-class postage paid at Lancaster, Pennsylvania. 


[Copyright 1959 by the University of Chicago] 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 21,1959 


A MATHEMATICAL INTERPRETATION OF MOTOR REACTIONS 
TO FIXED AND MOVING STIMULATION OF THE EYE 


MAXIMO VALENTINUZZI* 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO** 


The results of experiments performed by L. S. Stone and by R. W. 
Sperry on vision in salamanders are analyzed. A summary of Reichardt’s 
theory, which originated in the researches of B. Hassenstein on insects, 
is given, and an application to salamanders is carried out. 

Introducing a simpler mathematical development based on the one- 
factor theory of excitation, the same formal results are obtained. A com- 
plete model, which includes those results and tries to fit the salamander 
behavior described by Stone, is elaborated and translated into Boolean 
terms. As a complementary consideration, it is pointed out that the func- 
tional rigidity of the visuomotor system, deriving from a biological mo- 
lecular hysteresis, could be framed in the neuro-affinity principle, the 
spatial expression of which would turn out to be a biological field. 


Introduction. The investigations carried out by L. S. Stone and 
by R. W. Sperry have shown that it is possible to remove the eyes 
of the salamander and to transplant them either to the same animal 
or to another animal of the same or another species and that the 
eyes reacquire their vision after a certain time. A new retina is 
regenerated and a new optic nerve, which follows the pathway of 
the degenerated nerve, grows back through the optic chiasma and 
reaches the brain, where apparently it makes proper connections 
(Stone, 1947, 1955; Stone and Cole, 1943; Stone and Ellison, 1945; 
Sperry, 1943, 1945). 
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L. S. Stone (Stone and Cole, 1948; Stone, 1947) tested the func- 
tional quadrants of the retina in transplanted eyes under various 
conditions in animals having only one eye. He used either a colored 
bead impaled on the end of a white wire or a rotating black and white 
striped drum. The normal salamander moved toward and followed 
the bead as it approached all quadrants of the visual field. By 
means of the rotating drum, Stone verified that the normal animal 
moves in the same sense as the rotation only if it passes from the 
temporal to the nasal poles through the field of vision of one eye, 
and in such a way that this motion is counterclockwise for a right 
eye and clockwise for a left eye. 

Furthermore, Sperry (1943) and Stone (1947) discovered inde- 
pendently that when the grafted eye is rotated 180°, the salaman- 
der shows reversed responses. We refer the reader to the original 
papers for more details of the experiments. It is interesting to 
note that the reversal of the retinal quadrants appears at the time 
when retinal layers are beginning to differentiate (Stone, 1944). 

These researches prove that regeneration reestablishes the func- 
tional patterns of the retina, making corresponding connections in 
the central nervous system (Stone, 1947). One has also to think 
of a polarization within the optic tectum itself (Stone, 1944, 1947). 
In addition, transplantation to entirely different species suggests 
the concept that the eyes and the optic part of the brain are always 
similarly organized (Stone, 1953). In communication theory we 
would say that the vision system is a pre-addressed system. In 
general, these functional properties of salamanders have some 
common features with vision in the human (Stone, 1953). 

B. Hassenstein (1950, 1951) investigated the behavior of an in- 
sect, Chlorophanus viridis, when its eye was submitted to moving 
stimulation. Hassenstein employed, as a stimulating object, a ro- 
tating drum on the wall of which there were black and white stripes 
of a determined width and separation. By means of this device he 
could apply either a direct moving stimulation produced by the pat- 
tern or an indirect moving stimulation obtained by interference ef- 
fects using a second fixed cylinder with a system of slits. The 
animal reacted by following the sense either of the real or of the 
apparent motion. There was an optimal velocity of rotation which 
produced the strongest responses. For other interesting details of 
these researches we refer the reader to the original papers. 

Our purpose is to translate into a theoretical interpretation the 
physiological facts discovered by Sperry and by Stone. We believe 
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that the mathematical analysis of Hassenstein’s experiments de- 
veloped by W. Reichardt is particularly useful as a starting point. 

Hassenstein and Reichardt (1953) analyzed the sufficient and 
necessary conditions to be fulfilled by a transmission system in 
order to fit some functional characteristics of the central nervous 
system, and later (1956) they considered several possible struc- 
tures. Reichardt (1957) developed his analysis using a very sim- 
ple topological model. 

fetchardt’s theory. The model considered by Reichardt consists 
of two parallel lines of communication. It is structurally similar 
to the model of Figure 1 to be considered below. However, the 
model of Figure 1 lacks an element in each line and the properties 
of its elements are different from those of this model. At the af- 
ferent end of each line of the model under consideration there is an 
element (A or B) which corresponds to the transmitter of a com- 
munication system, and at the efferent end there is an element (S, 
or S,) representing the receiver. The subscripts A and B will be 
used to differentiate corresponding elements on the two equivalent 
lines starting with A and B respectively. The source of informa- 
tion is the visual field, and the destination of the messages is the 
musculature. For the beetle studied by Hassenstein and Reichardt, 
A and B represent two ommatidia. Each transmission line has 
three intercalated elements: D, or D, (differentiation elements); 
E,, E, (inertial elements); M ,, M, (multiplication elements), The 
lines are interconnected, there being a connection from D, to Mp 
and one from D, to M,. These lines constitute noise channels, 
and between A, B and S,, S, we have a black box within which 
the behavior of the actual optic neuron net should be solved. 

The luminance of the part of the stimulating object under ob- 
servation at time ¢ (e.g., a stripe of the drum) is expressed by L/(z), 
which plays the role of a stochastic variable with a probability 
distribution function p(L) and covers a range of values between 
zero and a maximum. If the stimulus object is moved in front of 
the elements A and B, these elements will register irregular changes 
L,, Lp in light intensity. These patterns traverse D,, Dp; trans- 
form into L%, L#; and p(L) becomes a symmetrical distribution. 
The message from D, reaches E, and Maz, and that from Dp goes 
toE, andM,. E, gives L4,,, and EF, gives Li... 

For convenience we may assume that the L* patterns are made 
up of unit pulses, represented by a 5-function. A weight function 
W(t’) is introduced, which is 0 for ¢ < 0 and a decreasing expo- 
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nential for “% > 0. With linearity of response and with L*(¢) as an 
arbitrary function, the expression 


epee i! ” W(L)L M(t - t’) dt’ (1) 


t 


holds. 

Now the calculation of the autocorrelation function of L* is 
needed. Considering the activity of the model at two different in- 
stants, such a function turns out to be 


, gO. oa 
PREM EC see ee @) 
where A¢ is the time elapsed between the stimulation of each 
one of the ommatidia, T, is a time constant, and 2, denotes the 
constant power spectrum. Expression (2) is obtained by the multi- 
plication L*(¢, )L,*(4,), by time averaging, and by application of a 
relationship between the power spectrum 1) (z#) and the autocorrela- 
tion function (Reichardt, 1957). 

Expression (2) says that between L4,(¢) and L#,,(¢+ AZ), there 
is a positive correlation which is symbolically written in the fol- 
lowing way: 

[L4,,.(¢), L4,(¢ + Az] > 0. (3) 


Considering that L#,, is functionally dependent on L*(¢), for a 
motion of a pattern from A to B we obtain the following from equa- 
tion (3): 


[L4,,(é), L#(2)] > 0; (4) 
and for a motion from B to A, we have 
[L4(t), L*,(2)] > 0. (5) 


At M, there is no correlation between the patterns L%,, and Ls. 

Reichardt (1957) analyzed the influence of the velocity of the 
motion. We briefly outline the argument. Consider a sequence A 
to B. The final outputs from $, and S, can be written as mean 
values of the products L4L%,, and LAL*,, respectively: 


THE [ff Lag, Ly LG. dd galy,s (8) 
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Here p(L 4, L}#..) is the joint probability of finding a value of L* 
in the interval (L%, L*%,+dL%*) and a value of L*,, in the interval 
(Lé.. UR. + ALR,,). The quantity p(L*, L*,,) has a similar mean- 
ing. Let us consider the element S,. Because there is no correlation 
between L% and L3},., the probability P(LA, L#,,) can be written 
as a product of individual probabilities. Furthermore, because the 
individual averages La and L#,, vanish due to the symmetry in 
the L * distributions, we obtain 


/LALR,, = La Le a (8) 


Let us consider the element S,. Since Lg and L*,, are cor- 
related, the joint probability p(L#, L%,,) cannot be written as a 
product of individual probabilities. To solve this case, the aver- 
age over a Statistical ensemble is replaced by a time average 
(Reichardt, 1957): 


oh +T/2 
L3L4,,> lim = L*(t)L* , (t)dt. (9) 
To3 -T/2 


When a signal L3(t) enters M, at an arbitrary we t, it becomes 
multiplied by the signal L¥,(t). This signal is integrated by a 
superposition of contributions which corresponds to equation (is) 
therefore we can write 


Eeo= [ W(t — tL *(t)dt” 2710) 


The only relevant value of ¢’ in the development of the super- 
position process is ¢— At, because ‘*At seconds ago,’’ the excita- 
tion represented by the symbol L¥(t) acted on E’,, leaving in this 
element a residual excitation which then would decay exponentially. 
Taking into account this fact and using equation (10), expression 
(9), with 

L4(t)=L4(¢—- At), (11) 
gives 
ot Ae 
= o*W(At) =—— e TE, (12) 
ug 


*7] * 
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where oc? is the quadratic dispersion or variance of the message 
L"(¢). 

Now we are in a condition to introduce the velocity of the light 
source. The time taken by one of the signals to move from A to B 
is At. Let us call Ad the angle described by the signal in the 
visual field. Hence we have 


abel (13) 


To approach more closely the actual biological system, Reichardt 
(1957) considers the D-element (with time constant Tp). The cal- 
culation of a time average like that of expression (9), as well as 
Shannon’s sampling theorem, gives 


2 LI 
DRT Ge (14) 
TUDE 


where & = T,/Tp. 
If tT, is introduced into equation (13), and the signals leaving 
S, and S, are combined, the total effect is given by the difference 


between the output of S, and the output of S,: 


2 mike o 
L§L¥,, - LAL§,, - —— . Bek *). (15) 


tr ea 
TUnTp 


An application to the behavior of the salamander. We now try an 
application of Reichardt’s theory to the motor reactions described 
by Stone. In the case of salamanders, the elements A and B could 
represent either two photo-receptors or two fields of the retina 
(right and left, upper and lower, respectively) of one eye. S, and 
S, would correspond to elements in the motor brain cortex with ap- 
propriate connections to some groups of muscles which are able to 
give deviations or movements of the animal toward the right side 
or toward the left side. 

We shall consider the reactions to fixed stimulation and to mov- 
ing stimulation. First we will analyze the case in which there is 
moving stimulation. 

The salamander with only one eye reacts when the motion of the 
stimulus is temporo-nasal. To take this experimental fact into ac- 
count, we assume that the model is asymmetric. Hence we modify 
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it by assuming only one interconnection; that is, L% to M, for one 
eye, and L 4 to M, for the other eye. In this way, the system should 
give a reaction to the sequence A to B, but not to B to A in the 
first case, and to the sequence B to A, but not to A to B in the 
second case. Obviously, correlations like those of expressions 
(4) and (5) hold, but now there is no functional symmetry. 

Let us next analyze the motor reactions to fixed stimulation. 
When an eye in the salamander is stimulated from the right side 
the reaction is toward this side, and when the stimulation comes 
from the left side, the animal reacts toward the left side. The re- 
actions follow a similar pattern for the upper and the lower parts 
of the visual field, and the mechanisms should be the same. 

We assume again the symmetrical topological models with the 
interconnections L% to M,, and L% to Mz, which, for the sake of 
simplicity, will be designated by (a) and (0), respectively, but 
now both of them are applied to the same eye (either right or left): 
(a) corresponds to the right retinal field, and (5) corresponds to 
the left retinal field. 

The stimulus acts at the same time on the photo-receptors A and 
B,. Signals L, and L, are transmitted by these elements and the 
following conditions hold: 


L 4(t) = L g(t); Ld) = LR); LA, (t) = 25,0). (16) 


For (a) a multiplication of the signals L4,, and L# takes place 
at M4; for (6) a multiplication occurs at M, of L4 and L#,,. In 
accordance with Reichardt’s analysis, the final outputs given by 
S, and S, can be written as mean values of the products roe Ee 
and L4L%_, respectively, and the expressions (6) and (7) are valid. 

There is no correlation between L* and L#,, at Mp for (a), and 
because of the same argument used to obtain equation (8), the out- 
put of S, reports zero excitation, and no motion toward the left 
side occurs. On the other hand, L*,, and L} are correlated, and 
now expression (9) applies. A positive output at S, should result 
and, in consequence, a motor reaction toward the right occurs. 

The signal L%,, is integrated by elementary contributions ac- 
cording to expression (10), but only for one value; that is, only for 
t—¢’=t*, the instant at which the stimulus acts, is there a con- 


tribution to the composition process. Therefore, 


L*, (¢#) = n W(t tL A(t)at’ = WEMEA(L- t%). (17) 
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Averaging over time, and considering equations (16), we obtain: 


DS 

i 

wx 
I 


+7T/2 
lim 1/T | L*(t)L*(t — t*) W(t*)dt 
To =e 


+T/2 
iia vr f L*2(t— t*) W(e*)dt 
-T/2 


T 00 


+T/2 
= W(t*) lim /t | L4*?(t— t*)dt 
T>0 
-—-T/2 


ee 
= o7W(t*) = — e 
TE 
For (b) we would deduce, by the same considerations, the fol- 


2. (18) 


lowing expression: 


LaLSs 


4TR,,=—e 7 >0. (19) 


Te. 

Expressions (18) and (19) mean that stimulation of the right reti- 
nal field should result in a motion of the animal toward its left 
side, and that stimulation of the left retinal field should give a 
motor reaction toward the right side. 


Application of the one-factor theory. The one-factor theory 
(Rashevsky, 1948) permits a description of Stone’s experiments as 
well as those of Hassenstein in terms of a very much simpler math- 
ematics, and under special conditions the final expressions are 
formally the same. We follow a suggestion by H. D. Landahl (per- 
sonal communication) and translate Reichardt’s topological model 
into a neuron model (Fig. 1). Here A, and B, represent photo-re- 
ceptors; A,, B,, A,, B,, A,, B, represent neurons. To simplify 
matters we shall suppose that A, and B, are of the ‘‘on”’ type of 
photo-receptor such that their rate constants a and } are fast com- 
pared with those of A,, A,, B,, and B, (cf. Householder and Lan- 
dahl, 1945, Chap. I). When they are suddenly stimulated by an in- 
tensity S, they will give a short burst of activity depending upon 
S. We shall suppose this dependence is linear, so that we may 
consider the burst of activity in A, as having an average intensity 
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MOVEMENT 
TO LEFT 


2 A3 MOVEMENT 
TO RIGHT 


FIGURE 1 


S4,, proportional to S,, and of a short duration ¢ The neuron Ae 
will be assumed to be of the simple excitatory type governed by 
the relation (Rashevsky, 1948; Householder and Landahl, 1945). 
dE, 2 
dhic 


APA, - €Ag), $4, = BSA, = BB,S,, (20) 


where f and f, are constants, and S4, is a measure of the ac- 
tivity of A, due to the stimulus S, acting on the A-chain. For 
convenience we neglect the thresholds of A,, A,, A,, B,, B,, and 
B,. We suppose that B, has the same parameters as A,. The ele- 
ments A, and B, are also assumed to be identical and to be of the 
simple excitatory type, having a’ instead of a for the rate constant. 
As in the previous case, we suppose that S, and S, are the same 
function of time, except for the displacement of A¢ in the time. We 
take the time to be zero when the stimulus first acts on the A-chain. 

Since S4,, and hence ¢4,, is appreciable only during the short 
time ¢, we find from equations (20): 


E4 = a6A,t=aBB,St, t<t 


app Sie *"~* ; = Gp Sen; t 2 ui (21) 


EA» 


where the subscript A has been dropped from S since the magni- 
tude of S, is equal to that of S,. Note that we have let ¢approach 
zero for convenience and that S’ is defined as 6B St by equation (21). 

The expression for €4, is obtained from an expression similar 
to equations (20), but with a replaced by a’, and $a, replaced by 
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PAg = Ba faze Since €4, = 0 for t= 0, if there has been no pre- 
vious stimulation, we then find, on integrating, 


é4, = aa’B,B,S’ (e-e "(a - a). (22) 

Expressions for €g, and €g, can be obtained in a similar man- 
ner, the only difference arising because of delay time At. 

The neurons A, and B, are assumed to have threshold values of 
two, so that only simultaneous firings of A, and B, will fire A,, 
for example. The input frequency from A, is proportional to €4,, 
while that from B, is proportional to €4,. The output frequency 
R, of A, is then, to a first approximation, proportional to the prod- 
uct £4, €B,; that is, if B, is the proportionality constant, we 
have, for ¢ > AZ, 


Ry stad BR, Bg eo eee le ee (a eae 


Fk, being zero for-¢ < Ad. 

Subtracting the similar expression for f,, the output of B,, we 
obtain &,,, the net output from the A-chain when the stimulus 
moves in the direction from A to B. We may suppose that the dif- 
ference is obtained by cross-inhibitory mechanisms (Landahl, 
1938). We then have (with ¢ > AZ) 


24° 2 02 
aap, BBs Ss “At At, ,-at—a’t 
ae (e? -— et ye at—a 


Rape (24) 
for the net output as a function of the time ¢>¢, R,, being zero 
for ¢< At. By integrating R,, over time, we may obtain the total 
excess R,, of R, over Rp. Integrating equation (24) from Az to 


co we find this to be given by the following expression: 
a? 2 42 
7 ee a vB, BB, S (en eA _ e7’Aty, (25) 


A 


If a=1/t,, @ =1/t), and a= Te/Tp, then it can easily be 
seen that equation (25) is the same as equation (15), except for 
the units. Thus under the above conditions, the neural network of 


Figure 1 yields a result similar to that given above by expression 
(15). 


A complete model for the salamander behavior described by Stone. 
In this section we give a complete neuron net model trying to fit 
the kind of stimulation and of motor reactions described by Stone. 
It implies fundamentally the type of communication lines consid- 


MOTOR REACTIONS TO STIMULATION OF EYE 337 


ered in the asymmetrical systems (a) and (6), but some intercon- 
necting fibers are inhibitory and, of course, the factor ¢ should be 
replaced by j in the respective mathematical functions. 

Figure 2 represents the model where the arrows S, and F,, re- 
spectively, summarize the inputs and the outputs obtained in the 
experiments. M, designates the receiver for motor reactions to 
the right side, and M, designates the receiver for motor reactions 
to the left side of the animal. 

52 
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FIGURE 2 


In accordance with the figure, we have the following table: 


Inputs Outputs 
a) Fixed stimulus from the right side: S; Ry 
Right b) Fixed stimulus from the left side: So Ro 
eye c) From right to left moving stimulus: S3 Rg 


d) From left to right moving stimulus: S3 0 
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a) Fixed stimulus from the right side: Sy R,4 
Left b) Fixed stimulus from the left side: Sx Rs 
eye c) From left to right moving stimulus: Se R¢ 

d) From right to left moving stimulus: S¢ 0 


The shifting of the stimulus on the retina is opposite to the 
actual shifting in the environment because of the geometric optics 
of the eye. 

For a Boolean description (McCulloch and Pitts, 1943; Rashevsky, 
1948) of this model, which seems to be useful, we refer the reader 
to the original report (Valentinuzzi, 1958). 


Discussion. The facts which have been referred to above in- 
volve several problems. 

According to Sperry (1948, 1945) there is an extreme rigidity in 
the functional organization of the visuomotor system of urodeles 
and this would seem to discredit an explanation of its functions in 
terms of relearning or other functionally adaptive processes. The 
reestablishment of reflex relations in the visual centers is quite 
independent of functional adaptation. These characteristics are 
implied in the word ‘‘polarization’’ used by Stone (1944, 1947). As 
we have pointed out above, the visuomotor system is a pre-addressed 
system. 

The physiological polarization of the retina is held fundamen- 
tally in the pigment cells, and the functional patterns in every 
sector of the retina become reestablished by cellular regeneration. 
Some intrinsic physico-chemical qualities in the ganglion cells of 
each quadrant of the retina may be assumed (Sperry, 1943). These 
specific qualities are permanently kept by the pigment cells near 
the ciliary apparatus (Stone, 1950a, 1950b). 

We may think of a sort of biological molecular hysteresis trans- 
mitted through mitosis. On this basic inherited hysteresis other 
remanent processes, such as learning, conditioned reflexes, etc. 
become superimposed. In this way the conceptions stated by N. 
Rashevsky (1948) and by W. S. McCulloch and W. Pitts (1943) re- 
ceive an enlargement. In reality, a living being turns out to be a 
maintained pattern (Wiener, 1954). From a pattern of this nature 
the functional rigidity of the visuomotor system derives. 

The reciprocal relations between the anatomical elements of the 
neuron structure of the vision system might be framed in a general 
concept, which could be called the neuro-affinity principle. 

The specific connections discovered by Stone and by Sperry, as 
well as those precisely described by Julia Apter (mapping of the 
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retina onto the colliculus, 1945), all of which occur either Sspon- 
taneously and naturally or after disturbing experiments (trans- 
plantation, rotation, etc.) should be related to some type of af- 
finity between the photo-receptors of special places in the retina 
and the optic neurons of special localization in the brain. The re- 
searches of P. Weiss give support to this line of argument (Weiss, 
1936; Shimbel, 1948). 

Stating the neuro-affinity principle, we would say: If two neurons 
possess a same feature they become functionally, and sometimes 
anatomically, related and can exchange information with each other. 
Impulse frequency could be such a feature. Weiss (1936) (Sperry, 
1943) considered a resonance principle. According to this con- 
cept, the basic functional organization of nervous centers is con- 
ceived to depend not upon a specificity of fixed neuron connections 
but upon the emission of qualitatively distinct excitatory agents 
and a selective receptivity to these discharges, permitting func- 
tional precision despite profuse and unrestricted interconnections. 

We may think of neuro-affinity as a special case of a more gen- 
eral principle which could be denominated cyto-affinity principle, 
and then we should say: Cells are related to each other if they 
possess some common characteristics. 

Finally, we wish to refer to the concept of fzeld implied by 
Sperry (1943) and other researchers. It is possible to imagine the 
situation of the cell relationships as phenomena deriving from a 
pattern which has the properties of a field. Every cell creates in 
its surroundings a biological field and, when the conditions of 
cyto-affinity are given, directing lines, such as lines of force, 
arise along which functional and, eventually, anatomical con- 
nections become established. 

From the cybernetic analysis carried on in this paper two sug- 
gestions emerge: (1) to construct an artificial experimental model 
with electric devices (automaton) which would reproduce the motor 
behavior of the salamander and prove or improve the assumptions 
introduced in the analysis; and (2) to design new quantitative ex- 
periments in salamanders. 


The author wishes to express his warm acknowledgement to Dr. 
Herbert D. Landahl and to Dr. Anatol Rapoport for their invaluable 
help. 
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This is a model for the time-variation of helium concentrations in lung 
wash-out curves. The helium (or other inert gas) is in a spirometer, 
which is connected by a common dead space to two separate dead 
spaces, each of which leads into a chamber. The chambers expand and 
contract, thus taking in some helium at each ‘‘breath.’’ Equations for the 
changes in helium concentration in each part of the system are set up; in 
this way difference equations are derived for the amount of helium in the 
spirometer after each breath, in and out, and complete solutions when the 
initial concentration is zero in all parts of the system except the spirom- 
eter. A simple solution when the chambers do not essentially differ 
(‘‘equal ventilation’’) is compared with the general case. The concept of 
“unequal lung ventilation’’ is discussed critically in relation to the 
model; some physiological interpretations are also mentioned. Numerical 
examples are given to show the effect of changes in various constants, 
in particular tidal volumes, end volumes, and the common dead space. 


1. Introduction. The degree of unequal ventilation of the lungs 
in emphysema patients, in asthmatics, and others is very important 
clinically and also patho-physiologically. One way of studying the 
process is to have the patient breathe a mixture of gases contained 
in a spirometer; this mixture consists of, say, 5 per cent of helium 
in oxygen or some other tracer gas not appreciably absorbed by the 
blood or tissues (Fowler e¢ al., 1952a; Becklake, 1952). The con- 
centrations of helium at successive breaths are measured. When 
the ventilation of both lungs is the same, this concentration de- 
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creases as an exponential function of the number of breaths, while 
when there are differences between the lungs, the rate of change 


is more complex. 


2. Models with and without dead spaces. This work can be 
looked upon as an extension of that of J. C. Gilson, P. Hugh- 
Jones, P. D. Oldham, and F. Meade (1955). They considered a 
two-chamber closed circuit system without dead spaces. The basic 
model is shown schematically in Figure 1A. B. F. Visser (1957) 
has extended this model by introducing two separate dead spaces, 
one for each chamber connecting it to the spirometer (Fig. 18). 

It is more realistic to assume that the separate dead spaces do 
not lead directly into the spirometer, but join to form a common 
dead space, through which air first passes on leaving and re-enter- 
ing the spirometer (Fig. 1C). This case has not been considered 
so far, probably because it is mathematically far more complex. 


FIGURE 1. Three two-chamber models: A with no dead spaces; B with 
two separate dead spaces; C with two separate dead spaces and a com- 
mon one. The largest and smallest volumes of the spirometer and the 
tidal and end volumes are also shown. 
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As will be shown, this common dead space can influence the de- 
gree of unequal ventilation considerably. Fortunately a different 
mathematical treatment from those used previously has brought the 
problem within manageable limits, and it is worth giving this in 
full. 

Readers interested mainly in the numerical examples and the 
physiological conclusions may omit this part. 


3. Details of the model. We consider Figure 1C. The two cham- 
bers have volumes F, and F, before inspiration, and as air is 
breathed in, they expand to F, +7, and F, + T,, respectively; 
corresponding to this, the volume of the spirometer contracts from 
H to H-~T, -—T,, so that the air is not compressed or rarified. 
The extra volume of gas reaching, say, the first chamber amounts 
to the tidal volume 7',; all the gas in the first dead space of vol- 
ume d, returns to the first chamber, as does a proportion, i,, of 
of the gas in the common dead space. 

On expiration, gas from the two chambers mixes in the common 
dead space and passes through to the spirometer, where it mixes 
uniformly with the rest of the gas in the spirometer. The idea is 
that in the common dead space the gas left behind after expiration 
does not mix with that in the spirometer at the end of inspiration; 
however, gas from the spirometer drives out all the gas in the dead 
spaces into the two chambers and fills the dead spaces, so that 
all dead spaces then have the same concentration as that in the 
spirometer. 


4, Notation. The notation looks arbitrary but was found by ex- 
perience to be convenient. 

The largest volume of the spirometer, at the end of expiration, 
is written as H. The smallest, equal to H-T, -T,, is written 
as equal to G - d, —d, ~d,, where: 


d, is the volume of the dead space leading to the first chamber; 


d, is the volume of the dead space leading to the second chamber; 
d, is the volume of the common dead space. 


Hence G is the volume of all these dead spaces plus the smallest 
volume of the spirometer (after inspiration). 

A lung is ‘‘well ventilated’’ at one breath if that breath brings 
in a volume of fresh gas that is an appreciable proportion of what 
is already there. Before the breath, this volume is F,, so a possi- 
ble measure of ventilation is T7/F. We shall see that the volume 
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of gas actually coming from the spirometer is more important. 
Clearly this is T, -d,-,d, = E,, say, for the first chamber. 
This corresponds te the effective Gaal volume, introduced by Gil- 
son et al. (1955, Appendix III). . 

The mathematical analysis pointed clearly to one particular 
measure of the degree of ventilation, namely to the ratio of EF, to 
the largest lung volume, F’, + T,. 

We write 

Tet Fels 
es ees 


E, 


Similarly, FE, = T,-—d, -—A,d, and 
d  eatheay 
2 E, 


The average degree of ventilation corresponds to 


d =5(¢; + do)- 


(A large ¢ corresponds, of course, to a low degree of ventilation.) 
Unequal ventilation of the two chambers is measured by 


“a 7 (d1 — $g)- 


The concentrations of wash-out gas, say helium, at the end of 
the r* expiration are given by 


u, in the spirometer, 
® in the first chamber, 
y, in the second chamber, 


2. in the common dead space, 


and the changes in concentration produced by the next breath in 
and out, are given by 


r 
X' aa, — &ry 
Y= p01 — Up» 
Z= 2,44 mt 


The proportions \, and A, must be assumed constant. If the 
lungs both begin and finish expending simultaneously, there can be 
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continuous flow only if 
Re etl eto) cand  \o mT oi(T, + To). 


This is not necessary in deriving the basic equations but will be 
assumed in the numerical examples. (This is justified in a model 
corresponding to ‘“‘parallel ventilation,’’? but not if there is ‘‘se- 
quential ventilation.’’ The evidence that the latter exists is in 
any case inconclusive. For references see Sec. 14.) 


5. The basic equations. The mass of helium entering the first 
chamber at the r + 1°" inspiration is given by 


@i,(F, +7,)—2,(F, +4,) =uF, + d,d,2,. (1) 
Similarly, for the second chamber, 
Vier (Fg +T.)-y,(F, + do) = UE, + 45A,2,. (2) 


The total mass of helium stays the same. We can write it as 
u,H; if at the beginning of the experiment all the helium is in the 
spirometer, u, is then its concentration. We have: 


uoH =2,(F,+4,)+y,(Fq +4.) +2,d, + u,H. (3) 
=241(F +71) + 9,41 (Fo +72) + 4,C. (4) 


Adding equations (1) and (2), it is easily seen that equation (4) 
is satisfied and so does not yield a new condition. 

During expiration, the mass of helium entering the spirometer is 
given by 
Hu,,,-(@-d,-d,-d,)u,=E,2,,,+ EoY,4,+(4,+4,+4,)u, (5) 
or, 


Hu, —-Gu,=E,2,,, + BoY p41) 


and at its end, the total mass of helium is 
UH = Hu,,, + (Fy +4,)t,4, + (Fo +4a)Yr41 + 4s2r41- (6) 


These equations are just enough to determine all concentrations 
uniquely, given initial values. 

This is not obvious, but can be seen as follows: if z,, y,, and u, 
are known, then equation (3) determines 2, Then equations (1) 
and (2) determine x ,, and y,,,, respectively; equation (5) gives 
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%., ame Ger equacer (@) ges ¢_.. Farter, all eqmacers @e 
cert equacmer 4.) Inve leer used. 
The cercem=scers im Ge stimemeter can be measured, = it ik: 
Test | ehmimake 2. ¥. ame 2_and cdeun 5 diiieence equation ime, 
i ts mack smpler lewewer & wort with Ge difkrerces (ere 
mares) X_. F_. Z_ ami U_. Thess ace sascy atauimed Oy ceaiscing 
> Oy ¢+ 1 im eae equader amé subaectiag. The equates Gen 


Ip is ateantageeas Gest @ sd equacems (2) am @) ang im eqae 
Gem (@) te replace + + 1 bye: She exqression Sea (F_ +€)+ 2, + 
@,) + 2,2, cam De Greeely substituted (Gime A, + A, = 2, elie 

Gu_==_.. @, + TL) +s... @, +7, -48% @ 
Simca, as is easily verified, © = #—E_ | 
We can new form four asic equacems: 
BE. Tatil, EE t)-2a- ee 
UF, + Fat, — 1, +2)- 247 -& , 
BG... — GU, +E 4, + EF. =& ee 
~o8_+ 622... + oF T= kh 
<a é 
& CQemamy She Gijfeeace equatium. We cam hem as mane 


i e 


4 CNEL SESTI_LatTAS 1F te we 
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The notation clearly shows the influence of unequal ventilation, 
however, the two coefficients of d,),r,/(E,£,) can be simplified 


since 
f(E, +E,)+ OE, -E,)=F,+T,+F, +T,. 


and = 
(¢ - 1)(E, +E,)+ OE, -E,) =F,+F,+4, +d,+d,. 
Further, 
¢? -@ = $,%93 
(fp — 1)? — 6? =(¢, ~- 1)(¢, - 1). 


This equation is symmetrical in that interchanging chambers 1 
and 2 makes no difference: @ then becomes —6, but A, and A, are 
also interchanged, as well as FE, and E,. 

The solution of such a difference equation is always the sum of 
two exponential functions of r and can be written 


U_=Aa’ + Bp’. (15) 
A and B depend on the initial conditions, and « and £ are the 


roots of a quadratic equation; this is found by replacing U,,, by 


A? and U, by A in equation (14). In fact the difference equation 
can be written as 


(A ~ o)(A - B)U,_, = 0, (16) 


where 
AU,_; + U,» A?U 


The values of Uy and U, determine A and B: 


U, =A+B; 
(17) 
Hence 
Ged 
a U,B _ BU, -U,. 
im ? 
i A B-a 
a B 
18 
paw l1— Wo oe 
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The concentrations are easily obtained from the differences: 


Up = Uy —(U, +U, +...4U,_,) (19) 
7 A(l-o") B(1- 8’) 
RG a ay eae 
A B Aa’ ; 
fh ii pe he ae Ty 


ok 1 Pad ano 1eaig 


After infinitely many breaths, the helium of course becomes uni- 
formly distributed throughout the whole system, so that its concen- 
tration is given by 

2! Huy 
See Mion MO did itd. itd.” 


(21) 


Clearly, too, « and 6 must both lie between 0 and 1; otherwise 
this equilibrium concentration would be infinite. It follows from 
equation (20) that 

A B 
L—4 - fp. 


(22) 


which provides a useful practical check on A, B, a and B, U,and U,. 

The concentrations in the chambers are best obtained in terms of 
the U, values already known. Equations (10) and (11) are solved 
as linear equations in X|,, and Y,,,. We have 


ys Hidalrs tex (ba 2A)GU; 


e+ 1 (23). 
E (dy — $2) 
es aA id taal (24) 
E'4(¢2 - $1) 
After many breaths, X, and Y, become, iin the limit, proportional 


to U.. 
To determine these limits, we choose a > A, then from equation 
(15) when r is large enough, U_,, ~ aU,. Hence 


Tr 


er Oe ae ; 
Tine Din = Lin 
Ui,, ree U, 100 U4 Ug, t Greg tees ree U, — t,, 


Hb, - {bq - 1)G/a3/{E (4, ~¢)}- (25) 
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Similarly, 


Lim a ine (Hd, - (6, - I)G/ahl/ {Ea ($2 - ¢,)$- (26) 


a) = 
re - ro u, Uso 


7. The initial conditions. The most likely initial conditions are 
Zy =VYo = 2% = 0; and we can put up = 1. Then from equations (1) 
and (2), 


a, =X, = 1/4; Y, = Yo =1W/dos (27) 
and from equation (5), 
Hu, -G =E,/¢, +E,/¢, =H(u, -1)+T,+T, -d, -d, -d, = 


H(u, -1)+E,+E,. 


=u, =Up = {6 (1-2) +8, (1 -2)} va (28) 
1 


From equation (3) we find that 


Hence, 


2, =Zy =h,/ gy + Ag/do- (29) 


For any other initial concentration u,(# 1), all these values have 
to be multiplied by u,. 

We must still obtain U, from the basic equations. First X, is 
obtained from equation (8), Y, from equation (9) and then U, from 
equation (10), after which Z, is obtained from equations (8), (9) 
and (11) with r put equal to 1. This will be shown in a numerical 
example. The general formulae are not very useful. 


8. Equal ventilation. We saw that the basic equations failed 
when ¢, =¢,; which corresponds to ‘‘equal ventilation.’’ This 
case is easily solved from equations (10) and (11) which now both 
involve FX 4, + E4Y,41,s Which can be eliminated, leaving 


AU,4, = GU ¢ be 1)/¢; 
so that 


-1)@)" 
U,= Uy {ea} . (30) 
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This relationship is in fact satisfied by equation (14) when @ is 


put equal to 0. 
But it is not the only solution; there is another, 


U, = Bp’, 


with different constants, 8 and B. The general solution is in fact 
still of the form of equation (15), even with 9 = 0. Clearly any ini- 
tial conditions consistent with the simultaneous equations must 
lead to B being zero. 


When uw, =1 and 2 =y,) = 2, =0, the general solution is very 
simple. Equations (8) and (9) become 


DX. (P= DX, + A,X. -—Z VE, = 0, 


with a similar equation for Y . 
Corresponding to (27), X) = Y) =Z,) =1/¢. Hence 


ype Eee 
ie Fe ae ay SR 
Mo 4 H 


(31) 
= pan 1-a-3, 
and X, = o/¢ = aX, = aYy; 
X =Y,=a7/¢; (32) 
X /U, =H/{(¢ - 1H - G)}; (33) 
2 =(1-a')/{d(1 - 3. (34) 


Numerical examples 
9. Two examples of equal ventilation: a. Fully symmetrical case. 
A reasonable set of values to take is: 


HAC Ces bi, om 1, dyn Bi 
G = 82, T, =16, F, = 44, d, =2, 


with every dimension in the second chamber equal to double those 


in the first. 
We have E, =2F, =12, 6,=¢,=5. It is easily found that, 
given u, =1, Ug = 0.144; Xo = Yo = 20 = 0.2; and that U, = 


354 M. E. WISE AND J. G. DEFARES 


0.144(0.656)’; further, X_ = Y, for allr, X, = 0.2(0.656)’, and 


1 Orie {1 — (0.656)' 
aa 1 — 0.656 
= 0.5814 + .4186(0.656)”, 


Similarly, #, = y, = 0.5814 {1 - (0.656)’$. 
6. Fully symmetrical case with larger dead spaces in relation to 
the volumes of the chambers. We can put 


H =100,d,=1,7, =4,F, =6, 
G =94,d, =2,T, =8,F, = 12, 
d, = 3, 


25 


and obtain: U, = 0.048; U, = 0.048(0.752)"5 ¥, = 37 + = (0.752)"5 


X, = 0.2 = Y,; X,= Y,= 0.2(0.752)% , = y, = 0.80645 {1 ~ (0.752)"}. 

The main difference is that the whole process is slowed up; 
even allowing for the higher end concentration, the rate of ap- 
proach to it is slower despite the higher ratio of tidal volume T to 
end volumes F. 


10. The exact condition for the concentration time law to be a 
single exponential. Using the example given below we shall study 
the condition forthe existence of a single exponential decrease of 
helium concentration with the number of breaths, i.e., ‘‘equal 
ventilation.”’ 

Example c1: Equal ventilation but the dimensions relating to the 
two chambers not in proportion. If we take F, =F, =50, T, = 20, 
PT, = 15, A, = 4/7; A, = 3/7; ds = 7; then ¢, = 70/(16 -d,); do = 
65/(12 —d,), and ¢, and ¢, are equal if 


65d, - 70d, = 200, 


which is certainly possible, e.g., d, = 4, d, = 6/7. 

In that case the changes in the spirometer and chamber concen- 
trations are still given by the simple formula. The results are 
shown in Table 1. ] 

This example clearly demonstrates that the relation Cs 
T,/F, is not sufficient to produce ‘‘equal ventilation.’”? In the 
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same way, the relation 7',/F, #T,/F, does not necessarily imply 
“‘unequal ventilation.’’ Only when the condition $1 = dq OF, in full, 


ewe = oars 
r, -d, — Ad, Ue -d, — ods 


is satisfied, is the difference equation of first order, which has 
always been taken to indicate equal ventilation. 

This important conclusion has not been formulated before. It is 
easily demonstrated that this condition can be extended to the 
case of n chambers in parallel. 


11. Complete numerical calculation for a case of very different 
ventilation ratios. Example d: We put H = 100, G = 93, r, =4, 
Se eg ae, fo 1, -d,-<1, -d,=1, d, =3,.A, = 1/3, dy = 
2/3, so thatE, =2, FE, =5, d, = 20, ¢d, =5, d = 12.5, 6 =7.5. 

Further, we put vu, =1, 2%) =¥) = 2) = 0. Then U, = 0.059, X, = 
0.05, Yo =0.2, Z=0.15, 40X, = 0.05 x 37+ 0.15 — 0.059 x 2 = 1.882. 
25Y, =2x18+2 x15 -0.059 x 5 = 3.605, whence X, = 0.04705, 
Y, = 0.1442, U, = 0.04672. 

In the difference equation the numerical values are given in the 
order in which they appear in equation (14). We have 


0 =100U,,, x 100 -U, {(193 x 87.5) + (7 x 2.5) + 5 - (62 x 65)} + 


980. . 116 — Gi x 58). 


r—1 


Obviously the two middle terms in the coefficient of U, contrib- 
ute very little, while the effect of the common dead space is small, 
but by no means negligible, in both the U, and U,_, terms. The 
most obvious effect of the common dead space again shows itself 
in the FE values and in changing the initial conditions. We obtain 


U.,, — 1.64767U, + 0.67084U,_, = 0. 
(A - a) (A- B)U,_, = 9; 
where « and £ are the values of 


0.8238 + 0.007864; «= 0.9125, B = 0.7352. 


The quantity inside the square root is small and comes from the 
difference between two nearly equal numbers. Hence these roots 
are rather sensitive to small changes in the coefficients. The nu- 
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merical values are: 


U 


0.059 

0.04672 
0.03740 
0.03028 
0.02481 
0.02056 


These values can be checked by obtaining U, from the basic 
equations. Z, is given by 


RPP wWDOHO|? 


3Z, = 100U, — 37X, — 18Y,. 
Next we find X, and Y,; we have 

40X, = 87X,+Z,-2U,; 

25Y, =18Y, +2Z, -—5U,, 

100U, = 93U, - 2X, - 5Ya; 


whence X, = 0.04398; Y, = 0.10343; U, = 0.03740, agreeing with 
the previous values. 
The general solution is 


U_ = A(0.9125)" + B(0.7352)". 


Using equation (18) we find: A = 0.01885, B = 0.04005; then using 
equation (20) the concentrations wu, in the spirometer are given by 


u, = 0.6330 + 0.2154(0.9125)" + 0.1516(0.7352)’. 
(Check: w,, = 100/158 = 0.63291.) 


The changes in concentration in the chambers are interesting. 
Using equations (23) and (24), we find that 


X 4, =(1000U,,, - 744U,)/60; 
Y 41 = (1767U, — 2000U,,,)/75. 


X,41 remains positive but tends to zero of course, as the U, do. 
The result for Y , the better ventilated chamber, is surprising, for 
it becomes negative after 18 breaths (see Fig. 2). In fact the he- 
lium concentration in this chamber overshoots its equilibrium 
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0-633! 


CONCENTRATIONS - as ratio of initial concentration 


° 5 ike) 1S 20 21 
NUMBER OF BREATHS r 


FIGURE 2. The concentrations at the end of the rth expiration for 
case d, in the two chambers (2, and y,) and the spirometer u,. All three 
concentrations approach that for uniform mixing, shown as a dotted line, 
asymptotically. 


(steady state) value and gradually returns to it. This phenomenon 
was first predicted by Visser (1957). 

Figure 2 shows the concentrations in the spirometer and both 
chambers compared with the equilibrium value. After 15 breaths 
the better ventilated chamber has only slightly less helium (.671) 
per unit volume than the spirometer (.689), while that in the other 
chamber is far less. 


12. How should the degree of unequal ventilation be assessed? 
In this and the next section, it is assumed that the two chambers 
do actually correspond to the two lungs. There is another less ob- 
vious interpretation too: ventilation in one lung and part of the 
other may be similar, but different from that in the rest of the 
second lung. 

The logarithm of the change in concentration between the r*> and 
r+ 1'*breath, U_, approaches the limiting straight line 


log U,~ log A +7 log a 
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and departs from it where the second component in the sum of two 
exponentials is appreciable. It is generally thought that the larger 
the departure from an observed straight line, the more abnormal the 
dimensions of the system compared with those for normal lungs. 
The problem then is how to measure and assess this abnormality. 
Experimentally it will be most conspicuous if B >> A and B is 
moderately different from 4 Then the early part of the curve will 
be dominated by the more rapidly decaying component. But if 
B << a, this component will decay after so few breaths that most 
of the log plot will still be linear. In any case the ratio B/A is 
more important than @/§ in indicating the degree of unequal ven- 
tilation. 

This degree can be interpreted in at least two other ways. In 
case d, B is larger than A and the build-up of concentrations in the 
two chambers is very different. Examples f and g are even more 
surprising. The concentration z, in the better ventilated chamber 
overshoots not only its equilibrium value (z,,, equal to u,.), but 
also, after this, overtakes and passes that in the spirometer! More 
precisely, it decreases more slowly than that in the spirometer 
(see Table 2). This clearly happens whenever the limit, as r — 
cc, of X /U_ is less than -1. 

The contrast between cases c2 and j is also interesting. In both 
cases the log U_y. r plot is almost wholly a straight line; in case 
j the concentrations z, and y, remain nearly equal to each other 
throughout, as they do (exactly) in cases a, 5, and c, but in case 
ce2, « differs considerably from y,. 


TABLE 2 
The concentrations for example /, at the end of the first few expirations 


in the better ventilated chamber (2,), the spirometer (u,), and the other 
chamber (Yr). 


r Lr ur Yr 

0 0 1 0 

1 0.509 0.591 0.257 
2 0.549 0.524 0.346 
3 0.535 0.506 0.394 
4 0.519 0.499 0.424 
5 0.508 0.494 0.444 
6 0.501 0.491 0.458 
tf 0.496 0.489 0.467 
oo 0.485437 
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A third approach is to calculate back from an observed curve to 
the physiological constants, provided of course that the curve can 
be reasonably fitted by a sum of two exponentials. This analysis, 
together with the actual concentrations (as well as their changes), 
would yield five constants, u,, 4, B, a, and £B, while H and T, + 
T, are also known. In a single experiment the small constants d, 
and d, might have to be guessed, and then there is more than 
enough information for estimating all the others. In a series of ex- 
periments of course everything could be tested. The results may 
be valuable even if the model is far from exact. 


13. Discussion of numerical examples in relation to experiments. 
In cases e to j the unit of volume could be about 30 cc. Ine andf 
all the basic variables are the same except d,. (@ of course is de- 
rived from the other constants.) By either the “‘log U_v.r” crite- 
rion or the ‘‘x, different from y,’’ criterion, the ventilation in case 
f is very unequal and also, of course, much more rapid than in e. 
Table 2 gives the first few concentrations. 

In recent experiments A. Bouhuys e¢ al. (1956, 1957) studied the 
effect of increasing common dead space and tidal volume together 
(but using an open system). In fifteen experiments a system with a 
common dead space of about 60 cc. and a total tidal volume of 620 
cc. was compared with a system with a common dead space of 
about 600 cc. and a tidal volume of 1100 cc. They concluded that 
ventilation became more uniform in the second case. 

From this point of view, cases e, f, and g can be compared. Of 
the three log U_v.r plots, that of e, which also has the largest 
common dead space, is certainly nearer to a straight line than g. 
In g, as in f, the two components are well separated, but in f the 
faster component disappears very quickly. 

The remaining examples show the effect of decreasing the end- 
expiratory lung volume F, + F,. Some authors (Mundt et al., 1941; 
Fowler, 1949; Cournand e¢ al., 1950) suggest that this improves 
uniformity of ventilation, while A. S. Gordon e¢ al. (1951) suggest 
the opposite. D. V. Bates et al. (1954) concluded from their ex- 
periments that it made very little difference. In our example A, all 
the constants are the same as in g except F, and F,, which are 
both 20 instead of 50. The build-up of x, compared with y_ is about 
equally different for A and g: in the log U , v- r plot for A the depar- 


ture from the straight line is well defined, but the approach to it is 
quicker than in g. 
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There is no doubt about the result in the extreme case j- Com- 
paring it with e, both have the larger common dead space but, for 
j, ¥, =F, =0, We see from the table that ventilation has become 
almost uniform in every respect. 


14. Physiological interpretation of the model. A more realistic 
model would consist of very many chambers instead of only two. 
About half of them, in one lung, might be connected to one sepa- 
rate dead space and the other half to the other one: these would 
still join to form a common dead space. There is little point in 
having three or four chambers instead of two; more realistic would 
be a statistical model of many small chambers in each lung, with 
different distributions of the physical dimensions, in particular T 
and F’, for each chamber and even different types of distribution 
corresponding to normal and abnormal lungs. 

Finally, a few words about conditions in the lungs that can lead 
to unequal ventilation. These conditions are frequently found dur- 
bronchospirometry and simultaneous separate determination of right 
and left F.R.C. (functional residual capacity) and may be due, 
e.g., to the presence of a partially obstructing bronchial tumor, 
unilateral increased modulus of elasticity, unilateral paresis of 
respiratory muscles (diaphragm), etc. 

It has often been stated that bronchial carcinoma is frequently 
accompanied by emphysema. The evidence is partly based on 
wash-out curves. From the above it should be clear that a second- 
order (or higher-order) semi-log plot may result from the unequal 
ventilation between the two lungs resulting from the bronchial 
tumor, so that, in this carcinoma case, the existence of emphysema 
should not be assumed on the basis of wash-out curves. 

The possible influence of sequential ventilation (Rauwerda, 
1946; Fowler, 1949) is obviously not accounted for in our model. 
The same applies to phase differences due to unequal time con- 
stants (Otis e¢ al., 1956). 

This work was begun at the Pneumoconiosis Research Unit in 
South Wales. One of us, J. G. D., would like to thank members of 
its staff for a very enjoyable and stimulating time there and the 
European Coal and Steel Community for sponsoring his visit to 
Britain. 
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Based on the author’s stochastic treatment of the unimolecular re- 
action which contains a mathematical rationale for considering random 
fluctuations inherent in the reaction mechanism (Bartholomay, 1958, 
Bull. Math. Biophysics, 20, 175-90), a new and simple formula is de- 
rived for calculating the rate constants of such reactions. This formula 
is then applied to data from the literature and the results compared with 
estimates obtained alternatively using current methods which attribute 
all irregularities about the usually expected smooth, decaying expo- 
nential curve to extraneous causes unrelated to reaction mechanism. 
Both estimates are seen to be almost identical despite the opposing as- 
sumptions about the origin of the random component of the kinetic data 
of such reactions. 


1. Introduction. In an earlier paper (Bartholomay, 1958) a mathe- 
matical model for describing a fluctuating time course for the uni- 
molecular reaction A—»B was derived, using the methods of 
probability theory. This stochastic model provides a feasible 
mathematical expression for that component of irregular fluctua- 
tions in the concentration-time data, which is traceable to the 
mechanism of reaction itself. In the present paper a new and 
simple statistical formula for calculating the rate constant, or 
equivalently the ‘‘basic probability parameter p,’’ based on the as- 
sumption of this model, is derived and applied to two sets of data 
on unimolecular reactions. 

The mathematical model which forms the basis of this paper is 
an extension of the classical empirical kinetic theory in the fol- 
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lowing sense. Whereas the earlier theory predicts unique concen- 
tration values n(¢) for A at any time ab initio, according to the 
formula 


n=nie ; (1) 


where n, =7(0), and K is the first order rate constant, in the 
stochastic approach this same formula has been demonstrated 
(Bartholomay, loc. cit.) to be the statistical mean value taken with 
respect to the probability distribution functions which define the 
stochastic model. Provision in the theory for uncontrollable ran- 
dom fluctuations about this smooth negative-exponential central 
tendency of the unimolecular reaction process forms the point of 
departure of the stochastic theory from the current deterministic 
theory. Because of the lack of a definite theory for analyzing 
fluctuations, the application of formula (1) to kinetic data by 
standard statistical methods has often led to ambiguities and re- 
sultant difficulties in comparisons of results obtained by different 
experimenters (see examples 1 and 2 in Sec. 3, for example). 

There would appear to be a need for a more extensive theory for 
the analysis of the macroscopic aspects of kinetic data. As a 
step in this direction, the stochastic model has been suggested as 
a means of isolating mathematically from concentration-time data a 
random component which is traceable to molecular events pre- 
ceding and accompanying the reaction, and therefore ‘‘inherently”’ 
associated with the unimolecular process. The existence of this 
source of random fluctuation had already been demonstrated and 
well substantiated in early work on extremely slow reactions such 
as radioactive decomposition (see Rutherford, Chadwick, and 
Ellis, 1951), where the range of variation is very large. The sto- 
chastic model suggests a quantitative criterion for measuring the 
extent of this variation in the general case. 

It is possible that such a model can provide an experimental 
basis for differentiating between extraneous and inherent fluctu- 
ations in data. Since this matter requires an extension of the 
original model, it will be discussed in a later paper. In the pres- 
ent applications, no attempt is made to partition the fluctuation in 
the data quoted. Instead, it is demonstrated by two examples that 
the assumption that the total randomness in the data is inherent in 
origin leads to values for rate constants which are well within the 
confidence limits for estimates of these same rate constants cal- 
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culated by standard statistical methods on the alternative assump- 
tion that it is due to experimental error alone. 

Section 2 contains a brief summary for easy reference of the 
pertinent equations of the stochastic model. Further details and 
the derivations will be found in the earlier paper (Bartholomay, 
1958). In Section 3, an appropriate statistical estimate of the 
basic rate parameter p is derived, using the ‘‘Method of Maximum 
Likelihood.’? The formula is then applied to the calculation of 
the first-order rate constant of some now classical data on the 
sucrose inversion reaction obtained by S. W. Pennycuick (1926) 
and to data of J. H. Raley, F. F. Rust, and W. E. Vaughan (1948) 
and on the gaseous decomposition of Di-é-butyl peroxide at two 
different temperatures. Comparisons with estimates by others and 
by other methods are also given. 


2. Principal features of the stochastic model for the unimolecu- 
lar reaction process. The mathematical model assumed in this 
work is a Markoff chain (see Doob, 1953) with continuous time 
parameter, with finitely many ‘‘states,’’ and with probabilities of 
transition from one state to another in a given time interval of 
length ¢ (the so-called ‘‘transition probabilities’’) which are in- 
dependent of the beginning point of that interval (the ‘‘homoge- 
neous’’ or ‘‘stationary’’ property). A ‘‘state’’ of the reaction 
system is defined by the concentration of reactant A expressed in 
integral units, the number 7 of molecules per given volume; i.e., 
where n, is initial concentration, ¢ seconds after the beginning, 
m can assume (n, + 1) possible values 0, 1,..., %), which define 
the states S,, S,,...5, . 

Certain mathematical assumptions were made which amount to 
the fact that if (n, ~ n) decompositions of A have occurred in the 
time interval (0, ¢), then in the subsequent infinitesimal interval 
(t, t+ At) the probability of a single decomposition equals [unAt + 
o(At)), where lim Ae = 0. On this basis the stationary transi- 

Aro At 
tion probabilities p,,(¢) giving the probabilities of transition 
S, — S, in time ¢ were derived (see Bartholomay, 1958, eqs. 


[7]-[21]): 


AC) a(t lero (ieh) (OSk<iSn,)), (2) 
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where ee is the symbol for the number of combinations of 7 ob- 


jects, taken & at a time. Thus, whatever the time of arrival of the 
concentration at state S,, this model permits the calculation of a 
definite probability of its arrival at a state S, in ¢ subsequent 
time units. It is this property, as will be seen in the next section, 
which makes the powerful ‘‘Maximum Likelihood Method’’ of sta- 
tistical estimation applicable to the estimation of the basic pa- 
rameter ,» (which corresponds to the ordinary first-order rate con- 
stant K shown in equation [1]). 

Alternatively, the question of what can be said in an a@ prior 
sense about predicting the future course of a reaction on the basis 
of a complete knowledge of only its initial state is also answered 
by this formula and allows a comparison with the classical de- 
terministic theory. Given the initial concentration n,, the proba- 
bility subsequently in time ¢ of any value n, where 0 <n <n,, may 
be expressed as: 


n 


p, (t) = ee) eilitet Guba Sodenthrs (3) 


Corresponding to this distribution, the ‘‘expected’’ or mean value 
of concentration n and the standard deviation, as functions of ¢, 
are easily calculated to be: 


E (n, t) = Moe (4) 


o (Nn, t) = Vn (e7#* — eH), (5) 


The identity of the right-hand side of equation (4) with the de- 
terministic function (1) (u in place of K) demonstrates that a sig- 
nificant consistency exists between the deterministic and the sto- 
chastic theories for the unimolecular reaction; viz., the predic- 
tions of the latter theory agree ‘‘on the average’’ with those of 
the former. The standard deviation (5) is therefore interpretable 
as a measure of the extent of inherent random fluctuations about 
values ‘‘expected’’ a priori, according to either theory. 


3. Derivation and application of a formula for estimating the 
stochastic rate constant ‘*y.’’ In the most general case, suppose 
that a unimolecular reaction has been repeated R times under 
identical conditions, and that in each run the value of the concen- 
tration of A has been measured at (S + 1) successive, equi-spaced 
times ¢, = 0, t, =, t, = 2,...,¢, =st,... t, =St=T over a 
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total interval of observation [0, 7]. Let Mm. be the value recorded 
in the 7th run at the time ¢, = st, where (r=1, 2,..., R); (s =0, 
1,..-, 8). The data therefore comprise an [RP x (S + 1)] ‘‘observa- 
tion matrix’’ 


m m , : oem : 


in which the elements in the rth row are the ordinates (concentra- 
tion values) of the rth sample curve produced by the process. If 
the reaction is unimolecular, the elements in each column will 
differ slightly from each other, but in random fashion, and the rows 
will each exhibit the same downward, negative exponential trend 
reading from left to right. Assuming the validity of the stochastic 
model described in Section 2, the ‘‘Likelihood Function’’ corre- 
sponding to M is given by 

L(m 


+> Mm, 


ID im = 


R Ss 
wv25 M5) = II II P (Ms, 5-19 Mp, 53 OR (7) 
1 


r™1 #«s 


where the multiplicative factor is the probability of a transition in 
time Tt from the concentration state m_ ._, to m, .. Formula (2) 


with m_ ._, replacing 7, m, . replacing &, and T replacing ¢ sup- 


plies this required information and leads, by substitution into 
equation (7), to the natural logarithmic expression: 


7 sal 


ois F log Rape ete 


ret ess 1 m 


(m,. et ten ms) log (-1 + ert \\.2 (8) 


The ‘‘maximum likelihood’’ estimate ji: of » is obtained by solving 
the equation 
d log L 
Op 


0, (9) 
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the result being the convenient arithmetic form: 


1 Rao Rove 5 
ae =p i rien] ~1 is phy Ms - (10) 


r=1 s=1 r=1 s=l 


This is the ‘‘general formula’’ for estimating the rate constant yp 
from R repeated runs each producing (S + 1) values, spaced T time 
units apart. In the special, common case when 4 is estimated from 
a single run, the matrix M has only one row and this formula re- 


duces to 
4. S Ss 
ii ae log a m._,| — log Da ™ ; (11) 
s=1l s=1 


where now the subscript 7 is omitted. The quantities (m.) may be 
expressed in any convenient units. 

For purposes of comparison with current methods of calculating 
the ordinary rate constant AK, an ad hoc assumption will be made 
that all randomness is inherent in origin, and formula (11) will 
then be applied to data on two different types of unimolecular re- 
actions which already appear in the literature. In each case, the 
substantiation of the unimolecularity property has received a great 
deal of attention. Results obtained using equation (11) will be 
compared with those using the earlier methods referred to in the 
introduction. In the first example, Pennycuick’s data (Pennycuick, 
loc. cit., Table II, p. 9) on the historically important sucrose in- 
version reaction could not be used directly because the readings 
were not made at equi-spaced times. The required equi-spaced 
values were approximated by curve fitting and interpolation and 
are given in the accompanying Table 1. All the results obtained in 
the present investigation are calculated from the data in this 
table. Pennycuick’s estimates and those of other workers which 
are quoted for further comparison are all based on the original 
data. In the second example, the original data (Raley et al., loc. 
cit., Table II, p. 91) on the decomposition of Di-t-butyl peroxide, 
a gaseous decomposition shown by Raley e¢ al. to be unimolecular, 
were used for calculations. 

Example 1: sucrose inversion. Applying formula (11) to the 
entries in the last column of Table 1 gives 


1 
fi or [log 1022.0 — log 931.5] = 3.729 x 107° min7?. (12) 
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DATA OBTAINED BY INTERPOLATION FROM PENNYCUICK’S 
TABLE II (¢ = 25 Min.) 


Time (Min.) % Initial Sucrose 

S e Remaining = Mm. 
0 0 100.0 
1 25 91.5 
2 50 83.4 
3 75 75.9 
4 100 69.2 
5 125 63.2 
6 150 57.5 
a 175 52.6 
8 200 47.0 
9 225 43.9 
10 250 39.8 
id 275 36.2 
12 300 32.8 
13 325 30.0 
14 350 PAP 
15 375 25.7 
16 400 2h 3k 
17 425 BAe 
18 450 19.4 
19 475 17.6 
20 500 15.7 
21 525 14.0 
22 550 12.8 
23 575 11.7 
D4 600 10.6 
25 625 9.5 
= = 1031.5 


If these same data are analyzed from the point of view of the 
earlier theory, using, say, the method of least squares, the ordi- 
nary rate constant K is estimated as a regression coefficient from 
the formula 


s 
Die (¢, — t) (log m, — log m) 
“ s=0 
K =—_________——__,, (13) 
Ss — 
eGo) 
s=0 
which becomes 
3411.06 a 
= 3.732 x 107°. (14) 


~ 914,062.50 
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It is interesting to note that the stochastic estimate ji is included 
in the 99 per cent confidence interval of K: 


99% confidence interval of K: [3.732 + 0.48] x 10: (15) 


Thus independent estimates of the rate constant based on two al- 
ternative approaches are almost identical, even though one is 
based on the assumption that random irregularities are an inherent 
feature-of every reaction process, the other, on the assumption 
that randomness enters only as part of the experimental procedure. 

The ambiguity of the latter approach is evident in the history of 
different treatments of Pennycuick’s data in which different meth- 
ods of data handling led investigators to varying conclusions 
about the unimolecularity of the sucrose inversion process, 
Pennycuick took his own data to mean that it might not be unimo- 
lecular and pointed to similar results obtained by several others 
who raised the same doubts. He carried out calculations which 
showed a drift in the K estimates calculated over successive in- 
tervals. The average rates of change for the seventeen successive 
intervals shown in his Table II covered a range from [3.613 x 10774 
to [8.745 x 107°], the grand average of these averages being 
[3.686 x 107°] min™*. 

G. Scatchard (1926), on the other hand, interpreted the drift in 
Kk values obtained by Pennycuick to mean that the estimated K 
was really a well-defined analytical function of time. Having de- 
rived an expression for concentration which involved an integral of 
k(t) as a function of time, he used it to correct Pennycuick’s 
measurements and concluded that the data were in as good agree- 
ment with the unimolecular theory as one could expect. 

L. J. Reed and E. J. Theriault (1931) re-opened this investiga- 
tion along purely statistical lines. They derived a new statistical 
procedure so extensive as to require, ideally, the use of a high- 
speed computer. The factors of trend, constant and varying ‘‘ac- 
cidental’’ errors were all taken into account. Applying this method 
to Pennycuick’s data, they concluded that his unusually precise 
data did, in fact, conform very closely to a reaction of the unimo- 
lecular type. They cautioned that the usual methods of calculation 
of rate constants may furnish misleading indications when the uni- 
molecularity property of a given reaction is being examined; that 
many of these methods are inconsistent with refinements in ana- 
lytical technique. For Pennycuick’s data of Table II, a rate con- 
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stant of [3.681 + .013] x 10°? min™!, which is almost identical 
with Pennycuick’s average estimate, was obtained. 

Example 2. The second example is the gaseous reaction, the 
decomposition of Di-¢-butyl peroxide, the stoichiometry of which 
has been shown to be (see Raley e¢ al., loc. cit.; Frost and Pear- 
son, 1953, pp. 30-31, 310-17): 


(CH,), COOC(CH,), — 2(CH,),CO +C,H,. (16) 


Making the usual assumption of the existence of a direct linear 
relationship between concentration and pressure produced by re- 
actant, the manometric data were given first-order kinetic treat- 
ment by Raley e¢ al. The data at two different temperatures were 
considered. Those obtained at 154.6°C. were used to calculate 
average rate constants over two sets of overlapping 3-minute in- 
tervals: 0~3, 2-5, 3-6, 5-8, .... A grand average for all of these 
was reported to be 


K = [3.22 + 0.08] x 1074 sec™!. 


154.6 


The data at 147.2° were used to obtain average rate constants 
over successive overlapping 20-minute intervals: 0~20, 2-22, 6-26, 
10-30, 14-34, ..., 26-46, the grand average being 


, —4 -1 
| aoe 19.45.2001 01.0.2. Se0G. -. 


A somewhat different transformation of the first set of data led 
A. A. Frost and R. G. Pearson (loc. cit.) to obtain a slightly 
higher estimate: 


¢ —4 -1 
K i546 = 18.46 + 0.07] x 10°" sec”. 


These experiments, also, have been much discussed in the litera- 
ture. However, the unimolecularity hypothesis and the original 
stoichiometric assumption seem to be tenable, it is generally 
agreed. 

Using just the values for 0, 3, 6, ..., 21 minutes at 154.6 °, the 
least squares estimate turns out to be K,,, , = 3.26 x10 * sec™'. 
And using 4-minute intervals on the 147 .2° data, the least squares 
estimate is K,,,. = 1.48 x 107* sec”'. The stochastic formula ap- 
plied to the data, partitioned in this way, yields the results: 


1 1000 olen 
s = — log —— = 3.28 x10 “ sec |; (17) 
Miss “789 © 943 
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1 1638.1 ‘ fi 
- els = 1.42 x10°* sec” (18) 
Miata ~ 949 °° 1588.4 


which are seen to be almost identical with the least squares esti- 
mates and with the values reported by Raley e¢ al. The stochastic 
formula applied to the final manometric data used by Frost and 
Pearson on 3-minute intervals gives the result [3.52 x 1044 sec}, 


which agrees closely with their result. 


4. Discussion. A new, statistically reliable, and very simple 
formula has been obtained for estimating the rate constant of a 
unimolecular reaction mathematically described according to the 
stochastic model introduced earlier by the author (1958). As 
equations (10) and (11) show, the application of this formula to 
concentration data, observed at successive, regular time intervals, 
involves simply summing the different concentration values, sub- 
tracting from this sum & the last term, taking the natural loga- 
rithm of the remainder; subtracting the first term from %, taking 
the logarithm of this remainder; and dividing the difference be- 
tween these logarithms by the common time interval between suc- 
cessive observations. 

If it is assumed that the main contributions to apparent random 
fluctuations about the central tendency of the data comes from the 
molecular mechanism of reaction itself, then this formula applies 
directly. If, on the other hand, experimental error is to be taken 
into consideration, further analysis is obviously necessary. How- 
ever, an interesting feature, stressed by the examples in Section 3 
is that in all cases these two sources of random fluctuations 
treated as alternatives, rather than additively, led to statistically 
indistinguishable estimates of the rate constants: as far as the 
value of the unimolecular rate constant was concerned, it made no 
difference whether the simpler formula (11) of the stochastic 
theory was applied, or other statistical methods and assumptions 
based on the deterministic theory, as exemplified by equation (1), _ 
were used. 

In Section 2 the point was made that the original stochastic 
model (Bartholomay, loc. cit.) predicts fluctuating concentration- 
time curves which agree ‘‘on the average’’ with the earlier de- 
terministic theory and which therefore incorporates the smooth 
predictions of the latter theory as a strong trend, its mean value. 
The close agreement obtained in these examples lends further 
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substantiation to the applicability of the stochastic model, es- 
pecially if it is noted that the rate constant formula (11) is a di- 
rect consequence of the main features of that model. This method 
of calculating the rate constant is ‘‘self-adjusting’’ in the sense 
that once a concentration value is observed at a given instant, the 
probability of obtaining the next value T seconds later is a func- 
tion of this present value, as well as T. Thus the formula, de- 
spite its apparent simplicity, is as fine-structured as currently 
employed successive averaging techniques. 


This work was supported by the Howard Hughes Medical Insti- 
tute. 
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In a previous paper (Rashevsky, 1959, Bull. Math. Biophysics, 21, 
299-308) we derived an approximate expression for the maximal speed of 
driving in terms of the reaction time of the driver. In the present paper 
the possible effects of unevennesses of the pavement, of such distracting 
stimuli as road signs etc. on the reaction time are studied theoretically, 
using previous developments of the mathematical biophysics of the cen- 
tral nervous system. In this manner expressions are derived which deter- 
mine the maximal safe speed in terms of road conditions and other dis- 
tracting stimuli. Effects of those conditions on fatigue are also discussed. 


In a preceding paper (Rashevsky, 1959a) we outlined a possible 
approach to the mathematico-biophysical aspects of driving. The 
latter brings forth a complex interplay of mechanical and human 
factors. As we remarked in the preceding paper, a number of stud- 
ies have been made on the effects of various physiological and 
psychological factors on the driving process. We mention here 
only a few, for example, T. VV. Forbes (1939, 1941, 1945, 1953); 
q. W. Case and R. G. Stewart (1956); H. J. Stack (1947). A great 
deal of work has been devoted to methods of testing drivers. Stud- 
ies of the behavior of drivers under simulated conditions have been 
made (Luce, 1958). None of those studies was, however, guided 
by a systematic theoretical approach to the biophysical aspects of 
the problem. We believe that here, as in any other science, the re- 
sults of experimentation will be much more fruitful when the exper- 
iment is guided by theoretical research which may indicate where 
to look for significant relations. It is not our purpose to develop 
here a theory of any given set of observed phenomena. Rather we 
propose to outline a general method of a mathematico-biophysical 
approach, which must of necessity precede any specific theory. 
Without adequate guidance by a theoretical picture we may find 
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situations where an apparently foolproof test is found to be mean- 
ingless (Buck, 1957). In the present paper we shall develop further 
some of the ideas suggested earlier (Rashevsky, Joc. ctt.). Famili- 
arity with the preceding paper is assumed. The notations used 
here are the same as those used before. 

We derived the following equation for the maximal possible ve- 
locity on an empty straight highway: 

ried erst er (1) 
OT 
where s is the width of the highway; s,, the width of the car; 6, 
the minimum safe distance of the car from the edge of the pave- 
ment; /,, the length of the car; 6, the average angle by which the 
direction of the car sometimes deviates from the true course; and 
tT, the reaction time of the driver. 

We shall consider first somewhat more in detail reaction time T. 
This is the reaction time to the stimulus which corresponds to the 
perception by the driver of a deviation of the direction of the car 
from the proper course. The reaction is the corresponding correc- 
tion made by a slight movement of the steering wheel. We shall 
regard, as a first crude approximation, the intensity of the stimulus 
as constant. In subsequent studies we shall drop-this restriction. 

If # denotes the intensity of excitation at the connection be- 
tween the afferent and efferent neuroelement (Rashevsky, 1948, 
chaps. XxXx-xxxii; 1959b, chaps. i-ili), if #, denotes the threshold 
of the efferent neuroelement, and 7 is the excitation time constant, 
then, neglecting the conduction time and considering purely excita- 
tory neuroelements, we have for reaction time T the expression: 


E 
At eh 2 
ae (2) 


1 
T = k log 
The expression is more complicated for mixed neuroelements, as 
has been shown by H. D. Landahl (Rashevsky, 1948; 1959b). Again 
we shall confine ourselves at present to the simpler case given by 
equation (2). 
For not too strong external stimuli S, the intensity E is given by 


E =B(S-A), (3) 


where £8 is a coefficient of proportionality and A is the peripheral 
threshold of the efferent neuroelement. If S >>, then E is simply 
proportional to S, We shall assume this situation here. 


AUTOMOBILE DRIVING 377 


In the mathematical biophysics of the central nervous system we 
frequently consider the situation where the neuroelements which 
mediate different stimuli mutually inhibit each other (Rashevsky, 
1948, chap. xxxi; 1959b, chap. iii). If m identical pathways which 
cross-inhibit each other are each stimulated with the same inten- 
sity, and if n is sufficiently large, the net result is an inhibition. 
This corresponds to our inability to react definitely to a situation 
in which too many stimuli of equal intensity are present. 

Any central inhibition is equivalent to an increase of the corre- 
sponding thresholds (Rashevsky, 1948, chaps. xxx, xxxi; 1959b, 
chap. i, ii). Hence if, in addition to the appropriate stimulus of 
the sight of a deviation from the straight course, there are other 
“‘distracting’’ stimuli present, the effect will be an increase of 
threshold & and hence an increase of reaction time T. We shall 
now study the effect of some important ‘‘distracting’’ stimuli in 
driving. 

First let us consider the effects of chocks and jolts due to un- 
evenness or roughness of the road. These effects are very much 
reduced by an appropriate design of the car, especially by proper 
suspension and sufficient weight. But on a sufficiently rough road 
the effects of roughness will be felt in any car. 

These effects are due to accelerations produced over different 
irregularities. They depend largely on the shape and size of those 
unevennesses. Yery roughly we may treat the effect in the follow- 
ing approximate manner. 

Let us consider the road surface as having a sinusoidal profile 
along its length, so that if 2 denotes the distance along the road 
from an arbitrary point, then the height y of the pavement varies as 


2 
y = acos — (4) 


where a and X are constants. The coefficient @ represents the 
‘‘height’’ of the ‘‘humps’’ which are spaced at a distance X. If 
a car travels along such a road with a velocity 7, we have @ = vb; 
therefore the vertical displacement of the car is proportional toa 
cos (27vt/X). The coefficient of proportionality o will be the 
smaller, the better the design of the car. Thus with 0<a<1, we 
have for the vertical displacement of the car as a function of time: 


Qrvt 
y = 4 a COS reo (5) 
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The maximum vertical acceleration at each ‘‘hump’’ is given by 


4 2v? 
ene oe 
or, putting 
40? ha 
ait tous (6) 
we have for the maximal acceleration 7, ., per hump, 
oe = vo”. (7) 


This acceleration acts as a ‘‘distracting’’ stimulus. It seems to 
be more natural to consider not the stimulation per hump, but rather 
the intensity of the continuous stimulation, expressed as the num- 
ber of individual short stimuli per unit time. Since per unit time 
the car passes over v/X humps, the number of those short stimuli 
per unit time is proportional to v. Each individual short stimulus 
is proportional to v? according to equation (7). Hence the inten- 
sity S, of the distracting stimulus due to unevenness of the pave- 
ment is proportional to v®. Denoting by c, the coefficient of pro- 
portionality we have 


Sy= oye" (8) 
From equations (6) and (7), and from the discussion above, it 
follows that c, is of the form: 


477 Xa 
ot = x3 be (9) 


Thus c, depends on the design of the car, through ; on the aver- 
age height a of the humps; and on the spacing X between them. 

Considering the case where S, is sufficiently larger than the 
perception threshold (a limitation which again will be dropped in 
subsequent studies), we find (Rashevsky, 1948; 1959a) that the 
central inhibitory effect of S, will be simply proportional to Sy. 
Hence the effect will be an increase of R by an amount proportional 
to c 1. Denoting by 8, another coefficient, we thus now have for 
the threshold of the efferent neuroelement, 


R=R,+8,c,v°. (10) 


Another stimulus, which may be assumed to vary as v3, is the 
intensity of the sound caused by wind. For this we put 


S, =c,v°. (11) 
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Assumption (11) is probably not correct, and is used here only 
for the sake of simplicity, as an example. The following arguments 
may be applied in the case of any other assumption. The final ex- 
pressions may merely become somewhat more complicated. 

The presence of signs along the road, or of other cars, which are 
spaced approximately uniformly, contributes to a distracting stimu- 
lus S,, which is proportional to » and may be put 


Nee Ba Us (12) 

For the time being we shall confine ourselves to the above three 
stimuli. Our purpose here is mainly to illustrate the mathematico- 
biophysical method of approach, and the above crude approxima- 


tions are sufficient for that purpose. 
Denoting by 6, and #8, two other constants, putting 


B.C, + B.C, =Ys 3°3 = (13) 


and considering the effect of all three types of stimuli, we now 
have instead of equation (11): 


R=R,+yv> +n. (14) 
Hence reaction time T now becomes 


: l i ae ae 15 

Rta ME rete ye, 7 Ge) 
Comparing equation (15) with equation (2) we see that the reac- 
tion time increases with velocity v. Therefore, as seen from equa- 
tion (1), the maximal safe speed is reduced by the roughness of the 
pavement and by the other distracting stimuli. If we substitute 
equation (15) into equation (1) and solve for v, we shall find the 
new expression for the maximal possible speed. A closed general 
solution cannot be obtained. A simple expression is, however, 
obtained for the case when E >>R, hy << yv> + nv, and where 
either the linear or the cubic term in equation (14) may be neglected. 


When E >> R, we can expand 


a E-R a ( =) 
auras og E air iea e 


and preserve only linear terms. This gives 


yu? + nv 


kE oe 
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If yu? >> nv, we have 


Vile 17 
i hae 


Introducing this into equation (1) and solving for v, we find, putting 
8-8, -26- 41, 
) 


4 
v= Vee (19) 
¥ 


If nv >> yv%, then a similar procedure gives 


AEk 
V= — (20) 
7 

As we stated, an increase in speed increases reaction time T 
and therefore reduces the maximal safe speed. If the driver desires 
to keep a given speed v which is over the safe limit determined by 
the above considerations, he can do it by increasing his attention 
to the driving. Neurobiophysically, this means the excitation of 
some appropriate center which results in an increase of EF, this in- 
crease compensating for the increase in threshold R. We may say 
that under the conditions discussed above, the driver is ‘‘rela- 
tively relaxed,’’ whereas if he must increase his concentration on 
the road in order to compensate for an increase in threshold R, he 
is ‘‘relatively tense.’’ We shall now derive an expression for the 
increase AE of attention, or concentration, necessary to drive at a 
given speed v. 

If for any given speed w the driving should remain safe, then 
equation (1) must be satisfied for that speed. Now, however, T in 
equation (1) is itself a function of the given speed v. Moreover, 
since instead of £, we now have E + AE, reaction time T is now 
given by: 


nA: (18) 


E+ AE 
= — log ——_________—.. 
k E+ AE ~(Ry + yv® + nv) 
Substituting equation (21) into equation (1), and using equation 
(18), we find, solving the resulting equation for AF: 


kA 
e° (Ry + yu? + nv) 
a PS AS Sal i 


‘ (21) 


AE = E. (22) 


e* sj 
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For v=o, AZ =o. Infinite speed requires infinite concentra- 
tion. For v = 0 we have AF = Ry -E. Introducing v = 0 and AE = 
Rk, -—E into equation (21), we find T =. In other words, for zero 
speed the reaction time may just as well be infinite, which is 
physically evident. The ‘‘attention’? AE increases with the 
speed v, 

The excitation AF cannot, however, increase indefinitely. It 
tends asymptotically to a saturation value AF, . (Rashevsky, 1948; 
1959b). The value AF,, is approached as » — . Since the right 
side of equation (22) increases indefinitely with v, beyond a cer- 
tain value v* of v, there is no value of AE which satisfies equa- 
tion (22). This value v* is the maximal safe speed for a ‘‘rela- 
tively tense’’ driver. This value v* is greater than the value given 
by equations (1) and (15). 

Hitherto we considered the instantaneous effects of distracting 
stimuli. There are, however, other effects which are of a cumula- 
tive nature. They all may perhaps be provisionally classed as 
“*fatigue’’ effects. Any sufficiently prolonged stimulus is likely to 
produce fatigue, which increases with time. Inasmuch as there are 
yet no ways of measuring the subjective feeling of fatigue, there is 
no point in deriving theoretical expressions for it. However, fa- 
tigue results in general in an increase of reaction times, which 
may be interpreted as an increase of R. We may therefore legiti- 
mately discuss possible cumulative effects of distracting stimuli 
on FR, and thus on T. The simplest way to do it is to assume that 
in addition to the instantaneous terms yo° and ny, the right side of 
equation (14) contains also terms which are proportional to time ¢, 
counted from the beginning of the driving, the coefficients of pro- 
portionality being themselves proportional to yv® and nv, corre- 
spondingly. Thus, with f and g as constants, we now put, instead 
of equation (14): 

R=R, + yv? + nv + fyv t + gnrt. (23) 


The ‘‘attention effort’? AE will in general also produce a fatigue 
effect pAEt, which should be added to the right side of equation 
(23). Since AE is a function u(v) of v, where u(v) is given by the 
right side of equation (22), we shall have in equation (23) a term 
pu(v)t. There may also be terms independent of v. Thus the dis- 
comfort of sitting in a relatively fixed position almost motionless 
at the wheel is likely to produce a fatigue effect, proportional to 
time. This effect will be of the form gt, where qg is a function of 
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the design of the seat. The more comfortable the seat, the smaller 
is g. Thus the final expression for F is: 
R=R,+ yo? + nv + fyv*t + gnvt + pulv)t + qt. (24) 
For the sake of simplicity we shall discuss here expression (24) 
only for a ‘‘relatively relaxed’? driver, for whom u(v) = 0. A simi- 
lar discussion holds for the general case. As ¢ increases, F will 
increase and will eventually exceed E+ AE,,.. At that moment T 
will-become infinite, and the driver will end with an accident. 
When R is near E + AE.., T varies very rapidly with R. Remember- 
ing that 


as? 


Py te (25) 


where L is the distance traveled, we see that the sum of all the 
cumulative fatigue terms in equation (24), for u(v) =0, is of the 
form 


Gi + a. fyo?) L. (26) 


Thus for a given distance traveled, L, the increase of RF as a 
result of fatigue depends on the speed v. If g = 0, or is very small, 
that is if the seat is very comfortable; if g is also small, that is 
the number of distracting stimuli along the road is small; then R 
increases as vw”. Since T increases rapidly with R, T will increase 
more rapidly than the square of the speed. For a given distance L, 
the slower we drive under those conditions, the better. If, how- 
ever, g and q are not negligible, then the expression in parentheses 


in (26) has a minimum for 
3 
1/ q 
VU = ad (27 
2fy 


There is in this case an optimal speed at which the fatigue ef- 
fects have the smallest value. As we see, the optimal speed in- 
creases with increasing discomfort of the seat. If the seat is very 
uncomfortable, the optimal speed may exceed the safe speed. In 
that case the driver may tend to drive too fast and thus take un- 
necessary risk. 

The problem of design of seats which cause a minimum of fa- 
tigue has been thoroughly studied in connection with fighter-plane 
seats (Hertzberg, 1955). Unfortunately, when sitting comfort is 
not related to the process of killing other people, the problem 
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seems to receive little or no attention. It may not be a bad idea 
to investigate the effects of sitting comfort on classroom perform- 
ances or efficiencies of conferences and the like. 

We shall now derive an approximate expression for the compo- 
nent 6, of 0, which is due to unevenness of the road. 

Consider the pavement as having per unit length nm) humps of 
height @ and width 6. We may roughly consider that every time the 
car hits such a hump, it will be thrown either to the left or to the 
right, by a distance proportional to ad. The coefficient of propor- 
tionality € will depend on the design of the car. It is a matter of 
chance whether the car is thrown each time to the left or to the 
right. Thus on the average, if there are no systematic asymmetries 
of the pavement, the car will be thrown with equal probability, %4, 
to either side. However, due to statistical fluctuations, it will oc- 
casionally happen that n’ consecutive hits all will throw the car 
in the same direction. The probability of this happening is (14), 


s 
If n’€ab equals the width ta 5 of the available pavement, then 


the car will come within distance & from the edge of the pavement 
(Rashevsky, 1959a), and the driver will have to correct the situa- 
tion by an appropriate turn of the steering wheel.* After a while 
the same situation recurs, and the effect is very much the same as 
the idealized effective zigzag line discussed in loc. cit. If n’ <<n, 
then the total number of possible consecutive sets of n’ humps is 
approximately n, per unit length. Hence the frequency of n’ con- 
secutive humps per unit length is approximately 


tN5 
2) Nos (28) 


But the car comes near the edge if 
8 
nab = 5-8, (29) 


*Actually the requirement that n’ consecutive hits throw the car in the 
same direction is too restrictive. There may occur occasionally more 
non-consecutive deviations in one direction than in the other, which will 
bring the car to the edge. Taking this into consideration complicates the 
problem considerably and will lead to different final expressions. We 
shall not discuss the exact problem here, inasmuch as the purpose of the 
following derivation is merely to show how the maximal safe speed v* 
may be expressed in terms of the unevenness of the road through 699. We 
have already seen how that unevenness affects v* through T. 
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or 
s — 20 


 9Eab- 


(30) 


Hence, the number of times which the car comes close to the 
edge of the road per unit length of the latter is 
Sao 
11 2£ab 
a=(5) sie 
Therefore, the average distance d traveled by the car between 
two approaches to the edge is d=1/m. In terms of the ‘‘equiva- 


lent’? zigzag line (Rashevsky, 1959a), the angle 6, (assuming as 
before that it is very small) is given by 


Tos (32) 


As a component of 6, the quantity 6, enters into equation (22) 
through equation (18). The coefficient € is a function of the de- 
sign of the car. To the extent that equation (22), as we discussed 
above, determines the maximal safe speed v*, we thus can express 
this speed in terms of the unevenness of the road both through 6, 
and through T. 


The author is indebted to Professor H. D. Landahl for a discus- 
sion of this paper and for useful critical comments. This work 
was aided by a grant from the Dr. Wallace C. and Clara A. Abbott 
Memorial Fund of The University of Chicago. 
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ON A LOGICAL PARADOX IMPLICIT IN THE NOTION OF A 
SELF-REPRODUCING AUTOMATON 
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The notion of automaton as used by J. von Neumann is formalized 
according to methods previously described (Rosen, 1958, Bull. Math. 
Biophysics 20, 245-60; 317-41). It is observed that a logical paradox 
arises when one attempts to describe the notion of self-reproducing 
automaton in this formalism. This paradox is discussed, together with 
some of the recent attempts to construct automata which exhibit self- 
reproduction. The relation of these results to biological problems is 
then investigated. 


The purpose of this note is to point out a paradox which arises 
when one attempts to provide a precise formulation of the notion of 
self-reproducing automaton, as originally introduced by J. von 
Neumann (1951) and discussed subsequently by a number of dif- 
ferent authors. If this paradox cannot be resolved, then the con- 
clusion to which we are forced is that the existence of a self- 
reproducing automaton is a logical impossibility. We shall in- 
vestigate the implications of this result for biological problems, 
and we shall also discuss in the light of this result several of the 
models for replicating structures which have recently been pro- 
posed. 

We must first make precise the basic notion of an automaton. 
As the term is used by von Neumann, an automaton is nothing else 
but what we have elsewhere (Rosen, 1958a) called a system (or 
biological system); henceforth we shall use these terms inter- 
changeably. In a subsequent investigation (Rosen, 1958b) we pro- 
posed a general method, based on the theory of categories, which 
is suitable for the representation of arbitrary systems (or automata). 
We shall suppose that the reader is familiar with this work. 
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According to our representation theory, every automaton may be 
regarded abstractly as a collection of sets and mappings which 
satisfy a number of simple conditions (cf. Rosen, 1958b, pp. 324- 
25). We are therefore justified in carrying through our subsequent 
discussion entirely within the framework of this formalism, re- 
placing the notion of automaton by that of a set of mappings, as 
described in loc. cit. In particular, one of the results of our analy- 
sis of this formalism (Rosen, 1958b, Theorem 4) showed that an 
arbitrary automaton could be constructed in a canonical fashion 
from single-output automata. From this result it follows easily 
that there is no loss of generality in restricting a discussion of 
the general properties of automata to the special case of the 
single-output automaton. For the sake of simplicity, therefore, we 
shall deal entirely with single-output automata in the sequel. 

As we showed (Rosen, 1958b), an arbitrary single-output autom- 
aton is to be represented by a mapping f:A — B, where A cor- 
responds to the set of admissible inputs to the automaton and B 
corresponds to the set of admissible outputs. The set A is in 
general expressible as a cartesian product of simpler sets (termed 
indecomposable in Rosen, 1959), which correspond to sets of ob- 
jects available directly from the environment of the automaton. 

We can now formulate in a precise manner the notion of a self- 
reproducing automaton on the basis of what we have said above. We 
have as yet made no restriction as to the nature of the elements of the 
set B which corresponds to the admissible outputs of the automaton 
represented by the mapping f above. If the elements of B should 
themselves be mappings, then the automaton represented by f is 
capable of synthesizing automata represented by the mappings of 
B. Itseems then that we can make the automaton f self-reproducing 
by requiring that f € B. 

It is precisely at this point that the paradox referred to earlier 
appears, for a mapping f is not specified until both its domain and 
range are stipulated. According to the definition of self-repro- 
ducing automaton given above, the range of the automaton cannot 
be stipulated until the automaton itself is given, since the autom- 
aton is a member of its own range. Thus, neither the mapping f 
nor its range can be specified until the other is given. In other 
words, the concept of self-reproducing automaton requires the 
concurrent existence of two pieces of information, each of which 
is logically anterior to the other. This paradox is apparently in- 
herent in any attempt to formulate the notion of self-reproducing 
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automaton; hence the very notion of such automata appears to 
contain an internal inconsistency. 

We must add a few words of qualification with regard to the above 
formulation of the notion of self-reproduction. As has been pointed 
out to the author by Harold White (personal communication), se! f- 
reproduction is usually recognized empirically by the replication of 
visible structural properties (i.e., in terms of what we have else- 
where (Rosen 1958a) termed fine structure). Our formulation, as 
proposed above, makes no reference to fine structure, but deals 
only with functional properties of components (i.e., only with 
what we defined in loc. cit. to be coarse structure). The above 
construction, therefore, is somewhat more inclusive than the 
customary one. For example, a given component /, which is re- 
alized by some definite fine structure, might conceivably produce 
as output a component which is represented by the same mapping 
f, but which possesses a different fine structure. Under our formu- 
lation, then, such a component would also be regarded as self- 
reproducing, although this situation seems to differ from the com- 
mon usage of the term. 

It may be instructive at this point to review von Neumann’s 
construction of a self-reproducing automaton and to observe the 
manifestation of the above paradox in his model. Von Neumann’s 
point of departure is an inference from a theorem of A. M. Turing 
(1937) regarding the existence of a universal automaton, which we 
shall denote by ®, which is capable of constructing an arbitrary 
automaton f when provided with a ‘‘suitable’’ description of f. 
In more formal terms, we have the following situation: let f: A — B 
represent an arbitrary automaton. There corresponds to each such 
automaton a ‘‘suitable’’ description of it, which we shall denote 
by Up. If S, is an admissible set of environmental materials, then 
the pairs (S, z,) are elements of the domain of ®, and we have 
®(S, 7%) = f. From this formulation it is quite clear that ® repre- 
sents a legitimate single-output automaton in the sense we have 
defined above. 

Therefore, by the universality property, ® itself must have a ‘‘suit- 
able’’ description, i. Therefore, if any element in the domain of ® 
is of the form (S5, 7), we are immediately in the situation of the 
paradox enunciated above. However, it might seem that we can 
avoid the paradox by stipulating that ® cannot be defined on any 
pair of which 74 is an element. This of course would force a modi- 
fication of the notion of universal automaton, but at first glance it 
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appears that such a restriction is sufficient to render the notion ac- 
ceptable and carry out the remainder of the construction proposed 
by von Neumann. 

That this is not so is seen from the following considerations: 
since the automaton ®, as modified above, must exist as a physical 
object in the real world in order for von Neumann’s program to 
have meaning, we must require that ® can be fabricated out of 
appropriate elementary objects by a finite sequence of operations 
in consonance with the laws of physics. That is, ® can be con- 
structed by an automaton W from a set S of admissible inputs. 
Since WY is an automaton distinct from ®, the definition of ® re- 
quires that ®(S,, 7,) = ¥. Now let us consider the automaton 
which we shall denote formally as (® + YW). This automaton is 
defined on the cartesian product S x S,, x z,,and obviously will 
produce the automaton (® + V) from appropriate elements of this 
cartesian product. But here again, we have a clear situation to 
which the paradox we have enunciated applies. Hence it appears 
that the notion of universal automaton involves the same type of 
difficulty as does the notion of self-reproducing automaton. 

We shall now undertake, in the light of the above comments, a 
brief discussion of the ingenious replicating structures described 
in the papers of L. S. Penrose (1958) and H. Jacobson (1958). 
These models fall roughly into two types. The first type, termed 
semi-conservative by M. Delbruck and G. S. Stent (1957), is a sort 
of mechanical analogue of the replication structure of the DNA 
molecule suggested by the model of J. D. Watson and F. H. Crick 
(1953) (see Penrose, loc. cit., p. 64 et seg.); strong experimental 
confirmation of the existence of this type of replication has been 
presented by M. Meselson and F. W. Stahl (1958). In this situation 
it is clear that since the replicating structure is composite and 
must be decomposed in the process of replication, there is no 
question of its functioning as an automaton, for an automaton by 
definition is left essentially intact by any process that it happens 
to implement. 

The second type of ‘‘replicating structure’’ is of a somewhat 
more subtle nature, although it may be considered as a special 
case of the first type of process. It will be more convenient to 
proceed at first in an abstract manner. Let A, B be sets, and let 
A x B denote their cartesian product. For convenience, we shall 
denote an element of A x B by ab, where a € A and 5 €B. Intui- 
tively, if we regard ab as a construct formed from the elementary 
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pieces a and 4, then we can envisage a mapping (automaton) which, 
when provided with the construct ab and the elementary pieces 
a and 5, can juxtapose the elementary pieces to form a new con- 
struct ab. That is, we have a mapping F:[A x B x (A x B)| > 
(A x B) x (A x B) with the property that F(a, 6, ab) = (ab, ad). 

From the above formulation, it is clear that the only automaton 
which occurs in the process is the one represented by the mapping 
F, to which the structure which is being replicated serves as an 
input (together with the elementary pieces of which it is composed). 
The domain and range of F are perfectly well defined, so that 
again there is no question of encountering the paradox enunciated 
above. It appears that recognition of the fact that the construct 
ab is not the active agent in the replication has been obscured by 
the simulation of self-reproduction of ad which is inherent in this 
type of model. 

Thus in the models of Jacobson (loc. cit.) we find that the struc- 
tures a, 6 and ab correspond to locomotive cars moving along a 
specially designed track; the track itself corresponds to the map- 
ping F. Similarly, in the model of Penrose (loc. cit.) we find that 
a and 6 are represented by mechanical shapes which are capable 
of interlocking when agitated horizontally on an enclosed track; 
here again it is the agitated track which is the actual automaton. 
Thus the functioning of these models is seen to be quite inde- 
pendent of the paradox described earlier and does not alter the 
conclusions we have drawn. In each case, the structure to be 
duplicated is not possessed of an active role in the duplication 
process but is utilized, as far as the theory of automata is con- 
cerned, in a relatively passive manner. 

We have seen that in the semi-conservative type of model the 
occurrence of our paradox is precluded at the outset by the nature 
of the duplication process. Nevertheless, it may be instructive to 
regard this type of duplication as at least a formal mapping process 
or automaton. We can write for such an automaton F, analogous 
to the above discussion, 


F :[B x (Ax B)x A] — (Ax B) x (Ax B), 


subject to the property that F(b, ab, a) = (ab, ab), where the ele- 
ment a occurring in the first co-ordinate of the range of F and the 
b occurring in the second co-ordinate are derived from the original 
composite structure ad. In this formulation, it becomes clear that 
both forms of replication procedures which we have mentioned 
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are essentially the same from a purely formal point of view; they 
differ primarily in what we have elsewhere (Rosen, 1958a) termed 
‘fine structure.”’ 

We may now inquire as to the biological implications of the 
discussion presented above. The most natural application of 
these results is, of course, to the process of the replication of 
biological materials, such as DNA. If a biological material cannot 
replicate itself in the fashion of an automaton (i.e., by causing 
the organization of elementary materials in its immediate neighbor- 
hood into a replica of itself), as seems established by the exist- 
ence of our paradox, then we appear to be left with the type of 
replication process exemplified by the models of Penrose and 
Jacobson. As we have seen, in this situation there exists an 
‘auxiliary’? automaton which can organize appropriate environ- 
mental materials into a copy of a suitable given structure, either 
by destroying the given structure or by having the structure play 
a role analogous to a ‘‘catalyst.’’ In the former case, since there 
is no inherent danger of encountering our paradox, this auxiliary 
automaton may be of a purely formal character (that is, the system 
in question may contain an induced mapping which will exist on 
the basis of the abstract structure of the system but which will 
not correspond to an identifiable physical landmark in an actual 
realization of the system), as in the quantum-theoretic model pro- 
posed by N. Rashevsky (1959). In this situation there exists a set 
of environmental conditions, derivable in principle from the laws 
of physics, whereby a molecule of a substance A in the presence 
of simpler molecules will undergo a reorganization, as a result of 
which the original molecule is destroyed and two new molecules 
of A are produced. We can regard this procedure abstractly as a 
mapping process, but this type of mapping is purely formal and 
does not correspond to the existence of an automaton, as we cus- 
tomarily employ the term. This type of situation has the advantage 
that we do not need to investigate the question of whether and how 
the auxiliary automaton is replicated (since we deal in this case 
not with a physical structure but with a set of definite environ- 
mental conditions). 

Finally, it may be useful to examine the relation of the results 
enunciated above to the notion of ({l,®)-system, as introduced by 
us (Rosen, 1958a, 1958b, 1959). The (M,&)-system was primarily 
constructed to provide a model for metabolic activities of a single 
cell, but as we have remarked, it is necessary for such a theory 
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to deal with the replication problem if it is to have any claim to 
completeness. In one of these works (1959) we proposed a model 
for the duplication of the components F, which correspond to the 
genetic material of the cell, based on the notion of induced map- 
pings which exist formally in the abstract representation of the 
(N,8)-system. It will be noticed that these induced mappings 
play a role very similar to the auxiliary automata which appear to 
be required in order to fabricate structures which can be duplicated. 
It will further be seen that since the domains and ranges of all the 
induced mappings are perfectly well-defined sets, there is no 
question of encountering the paradox we have described above. 
Since these induced mappings will correspond to environmental 
objects (at least in the form in which they were enunciated in 
Rosen, 1959), there is no need to consider mechanisms for their 
replication, so that we avoid thereby an infinite regress which 
would destroy the usefulness of the model. 

Nevertheless, it appears that the replicative mappings induced 
in abstract ({,#)-systems in the manner described (Rosen, 1959) 
may not correspond to the actual coarse structure of biological 
systems. In order to conform to the discussion set forth above, it 
seems that the mappings which are of importance for replicative 
procedures should be chosen so as to be of the form of the auxiliary 
automata occurring in the models of Penrose and Jacobson. This 
modification does not appear to be too difficult, and appropriate 
investigations are currently in progress. 


This work has been aided by United States Public Health Service 
Grant RG-5181C2. 
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MATHEMATICAL BIOPHYSICS OF COLOR VISION: 
Ill. COLOR CONSTANCY 
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A mechanism is presented which can account for certain aspects of the 
phenomena of color constancy. The mechanism involves interaction be- 
tween a given region and the remaining field. Each region is represented 
by a color center having the structure previously introduced (Landahl, 
1952, Bull. Math. Biophysics, 14, 317-25) to account for a number of 
phenomena of color vision. The trichromatic, symmetric mechanism is in- 
troduced for simplicity. The interaction is such that collaterals from each 
of the primaries representing the background send elements to each of the 
centers corresponding to the primaries representing the spot. However, 
the collaterals impinging upon unlike centers are excitatory while the 
collaterals impinging on like centers, corresponding to the same primary 
colors, are inhibitory. With proper choice of coefficients, the result is 
that for small changes in illumination, the resulting apparent color is un- 
changed. However, for greater changes in the color of the illumination, 
there results a distortion of the apparent color. A number of examples 
are illustrated numerically. 


The observation that the colors of objects tend to remain rela- 
tively constant under rather marked changes in the color and in- 
tensity of illumination is very old. The effect can perhaps be most 
easily demonstrated by viewing a large field containing colored ob- 
jects through tinted glasses which are covered on the periphery so 
that all light must pass through the glasses. H. Helmholtz (1867) 
and E. Hering (1878) considered the approximate constancy of color 
to be one of the most significant and interesting facts of the psy- 
chology of vision. Helmholtz emphasized the role of experience in 
explaining the phenomenon. Hering, while pointing out the role of 
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such factors as retinal adaptation and the accommodation by the 
pupils, nevertheless introduced the concept of ‘‘memory color’’ in 
accounting for the observations. D. Katz (1911, 1930) stressed the 
importance of the general illumination of the visual field. 

Since some fish appear to show color constancy to an even greater 
extent than does man (Burkamp, 1923) and young children under 
many conditions perform about as adults (Burzlaff, 1931), one might 
take the view that the contribution due to learning may not be the 
most important component, an important part of the effect being due 
to interactions at a low psychological level. Although this is not 
essential for our purposes, such a view is consistent with the model 
to be presented below. 

Another related phenomenon is that of color contrast. If a test 
patch is placed against a larger uniformly colored background, the 
color of the test patch tends to move away from the color of the 
background. If the spot is gray, it may take on the hue of the color 
which is complementary to the color of the field. 

It is the purpose of this paper to consider a modification of a 
neural mechanism for color vision (Landahl, 1952, 1956) which 
shows the property of approximate color constancy and which can 
account for some aspects of color contrast. 

Consider first the situation in which a spot of light is observed 
against a uniform background covering the entire visual field. Let 
R(A) give the reflectance of the spot as a function of wave length 
A, #(A) being that of the background. Let /(A) = /,7(A) characterize 
the illumination. Let p,(d) be the sensitivity of tie i*® primary as 
a function of wave length. We shall use the linear, symmetric 
mechanism involving three primaries (Landahl, 1952). The spot 
— act as a stimulus S having three components S,, S,, Ss given 

¥ 


go 29) [ p,(A) R(X) t(A) da. (1) 


For the background the corresponding S, is given by equation (1) 
with R replaced by R. 

When natural white light is used, for which /(A) = 1, 2(A), then 
S,; becomes S$,’ and 8, becomes Ae For convenience, we shall sup- 
pose that the response, the total tate of firings of neural ele- 
ments, has the same numerical value as § so that the units of p are 
adjusted in such a way that 


S,=1,/3 when R =1 and (A) =1,i’()). 
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From the properties of the mechanism for color perception given 
previously (Landahl, 1952) a spot of color illuminated by white 
light, represented by S’, would result in excitation in the nei ghbor- 
hood of a given point in the color region C,C,C, of Figure 1 (also, 
ef. Fig. 1, Landahl, 1956). We shall denote this point by S’. If the 
illumination is characterized by /(A) other than white light, then 
we shall denote the point by $. The X, coordinate of this point 
will be given by [Landahl, 1952, eq. Be with y, = yj: 

y S,- 28, -*8, 

Lys e cg ©) 

i j k 
For convenience we shall also write 
S = X (8) (3) 

so that the point $ represents the color of the light coming from the 
spot and impinging on the eye while S represents either §(A) or the 
components S,, S,, S, of the light reflected from the spot. 

The color of any illumination /(\) can also be represented by a 
point in the region C,C,C,, it being at that point corresponding to 
the color produced when the spot illuminated has constant reflect- 
ance, provided that the background is white. We shall use the sym- 
bol I to represent the point corresponding to the illumination /()). 

We next consider the interaction between the spot and the sur- 
round. For simplicity we consider only the effect of the surround 
upon the spot and suppose that the interaction is symmetric with 
respect to the primaries. Let the center C, corresponding to the 
spot be inhibited in proportion to the intensity of excitation at the 
center Cc corresponding to the background, but let C, be excited by 
éléments from C,,% #j. The inhibitory elements from C, go only to 
C, while the excitatory elements from C: go to both C, and C;,.- 
Thus let the coefficients be 2/3 and 1/3 respectively. Let the ap- 
parent intensities of the spot and the background be equal, that is 
=S,=5S,. Then instead of an excitation S; at the center C,, 
there will be a net excitation o, at C,, given by 


o,=8,- 78; + 75; + 5, =S8,-8;+ 5 2S;. (4) 


The point of maximum excitation in C,C,C, will now be at a point 
which we shall label o, instead of at a point $. The X, coordinate 
of the point of maximum excitation will then have the value 


(S, a S,) rz aS; =< 8) i + (S; fs S,) (5) 
SS S, +8, +8, : 
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or also because of equations (3) and (4), 
oomik(S 8 40d sSeiale (6) 


Let the surround be white so that R = 1 and let the illumination 
be nearly white so that 


I(x) = 1, (A) + 1,8) (7) 


where & is uniformly small. Let the spot be fairly close to white so 
that 


R(A) =1-A() (8) 


where A is also uniformly small. Under white illumination the spot 
would have the component S,” given by 


a = [e,(ig# O)[L-AQNIAA. (9) 


t 
The point on C,C,C, corresponding to this color would be given by 
S’ = X(S’). Under the changed illumination the spot will be repre- 
sented by S = X(S) in the absence of interaction with the boundary. 
However, due to the interaction it will be represented by a point 


o =X(S-S+ = =S;,). (10) 

But since z’ (A) represents white light, expression (9) can be written 

S/=i1,- In [PEO AAA, (11) 

while the component o, of o from equation (10) becomes 
0,= 41, ~1, [ pQ)# (A) Q)aA 
—T, JP, ASA)AQ)AA. 


Since the last term of equation (12) involves the product of 5 and A, 
it follows from a comparison of (11) and (12) that the model gives 
color constancy to a first approximation. 

In order to illustrate the effect of larger changes in illumination 
we shall give the results of some numerical calculations. In the 
first example, shown in Figure 1A, a spot is observed under six dif- 
ferent illuminations. The spot, S’, is assumed to be an unsaturated 
primary when viewed under white light against a white background. 
In the upper line of Figure 1A the illumination (I) has the same 
color as that of the spot but is more saturated. The reflected light 
viewed against a white background would appear to be almost satu- 
rated (S) but because of the interaction with the background, the 


(12) 
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apparent color, a, is actually considerably less saturated than when 
observed under white light with a white background. This effect is 
observed with various colors. In the second line the result of mov- 
ing the illumination toward white (0) is shown. The apparent 
color moves towards S’. In the third line, in which the illumina- 
tion takes on a slight hue of the spot color, the color appears to be 
slightly more saturated. When the illumination is approximately 
white, as in the fourth line, g and S’ coincide. As the illumination 
is moved toward the complementary color of S’, the color appears to 
be less saturated. In the last line, where the illumination is fairly 
saturated so that the reflected light (S$) also has the hue of the 
complement, the color appears to have the original hue but is much 
less saturated. 

In Figure 1B the same six illuminations are applied to a less 
saturated color, the color having again the hue of a primary. The 
results are very similar to those of the previous case. 

In Figure 1C the same illuminations, with one additional, are 
used with a spot of color which is an unsaturated color having the 
hue of the complement of the primary. In this case, as in the previ- 
ous cases, the color always appears to have the original hue, but 
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an illumination which is too saturated in either extreme tends to 
reduce the saturation. But here a small amount of the primary added 
to white in the illumination acts to enhance the saturation of the 
spots which has the hue of the complement. 

In Figure 2 are shown the results of five different illuminations 
on each of two different colors. Subscripts are added so that the 
corresponding points /,, $, and a, can be located. The illumina- 
tions chosen are two unsaturated colors having the hue of a pri- 
mary, two unsaturated colors which are complementary to this pri- 
mary, and white. In the upper part of the figure (at D) the spot, 
viewed under white light, is an unsaturated primary different from 
that of |, The values of S fall along a line containing C, and S’. 
The apparent colors, however, move along a curve in which the 
saturation increases and the line changes in an opposite direction. 
In the lower part of the figure (at #), in which the spot is an un- 
saturated binary under white light, the result is similar. The spot 
appears more saturated and changes color as the illumination 
changes from a fairly saturated primary through white to a rather 
saturated binary. 

In the above model we have assumed specific numerical values 
for the coefficients of interaction. In this case when the entire 
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field has uniform reflectance and is illuminated uniformly by any 
color, no color should be observed according to this model. If the 
coefficients of interaction are all somewhat smaller than have been 
assumed, then illuminating the field with a strong color will result 
in the entire field appearing to be of a pale color in spite of the 
rather saturated color of the illumination used. This situation 
corresponds more nearly to reality. If the coefficients are slightly 
different for the various primaries, then some saturated colors will 
be more and some less pale when the entire field is uniformly il- 
luminated. It may be possible to obtain restrictions on the numeri- 
cal values of the coefficients of interactions if such observations 
can be made quantitatively. 

Thus far we have considered that the intensity of the background 
is the same as that of the spot being observed. In this case we 
have not considered the effect of changing the intensity of the spot 
or background. In the particular model as given, the effect of in- 
creasing the intensity of the surround is to cause the point repre- 
senting the color of the spot to move further away from the point 
representing the illumination. Since this does not occur to any ex- 
tent, it is more likely that it is the relative intensity which enters 
expression (5). If the logarithmic model (Landahl, 1952) were used, 
the effect of intensity would be expected to be less important. 
Consideration of this case would be of interest as would be the 
effect of the introduction of thresholds of the various elements. 

In the above case we have considered only a uniform background. 
In general, one would expect that the effective background would 
depend upon some sort of integral which in turn depends on the way 
in which the interaction is related to the distance between the spot 
and any selected region. If there are many colors present, the ef- 
fective background will depend upon the spatial distribution of in- 
tensities of the colors. The effective background for one region 
may be different from that of the next. This will lead to color con- 
trast effects. It is convenient to refer to color constancy as the 
effect due to the over-all interaction while referring to color con- 
trast as the effect arising from the contribution of the interaction 
at short distances, perhaps being primarily border effects. The 
problem is then to determine the nature of the interaction from ob- 
servations on color constancy and color contrast. Such effects as 
the disproportionate role played by the brightest part of the field 
or the presence of a white area must be taken into account in the 
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interaction. The effects in which more than one color can be seen, 
depending upon the set (e.g., color memory) or upon the way in 
which the pattern is integrated, will require separate consideration. 
"lowever, if it is possible to obtain a function to represent the ef- 
fective background for a given spot, then by using a modification 
of the above mechanism it may be possible to account for a number 
of phenomena of color constancy and of color contrast with a single 


model. 
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NOTICE 


1960 INTERNATIONAL CONGRESS FOR LOGIC, 
METHODOLOGY AND PHILOSOPHY OF SCIENCE 


An International Congress for Logic, Methodology and Philoso- 
phy of Science will be held at Stanford University, Stanford, Cali- 
fornia, U.S.A., from August 24 to September 2, 1960, under the 
auspices of the International Union for History and Philosophy of 
Science. 

The proceedings of the Congress will be organized into the fol- 
lowing eleven sections: 

1. Mathematical logic. 

2. Foundations of mathematical theories. 

3. Philosophy of logic and mathematics. 

4, General problems of methodology and philosophy of science. 

5. Foundations of probability and induction. 

6. Methodology and philosophy of physical sciences. 

7. Methodology and philosophy of biological and psychological 
sciences. 

8. Methodology and philospphy of social sciences. 

9. Methodology and philosophy of linguistics. 

10. Methodology and philosophy of historical sciences. 

11. History of logic, methodology and philosophy of science. 

The proceedings will consist of a number of invited addresses, 
in addition to brief contributed papers. The closing date for sub- 
mission of abstracts of contributed papers is March 1, 1960. 

Information about membership fees and other details of the con- 
gress may be obtained by writing Professor Patrick Suppes, Serra 
House, Stanford University, Stanford, California, U.S.A. 
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